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flat  plate  to  a  typical  air  shock  wave  were  investigated. 
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material  and  geometric  nonlinearity  were  studied. 

In  this  study,  (1)  the  necessity  of  the  modification  of 
old  armored  vehicles  was  reviewed  and  (2)  the  NASTRAN  code 
was  employed  in  this  investigation,  (3)  the  theoretical 
background  of  the  nonlinear  transient  analysis  was  described, 

(4)  a  step  by  step  procedure  of  analyzing  the  dynamic  load 
problem  by  shock  wave  using  NASTRAN  code  was  developed,  and 

(5)  sensitivity  analyses  were  performed  and  also  difficulties 
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I .  INTRODUCTION 


Conventional  tanks  (old  tanks)  which  have  steel  armor, 
can  be  completely  defeated  by  one  round  of  modern  projectiles 
such  as  anti-tank  missiles.  That  has  been  proved  in  the  1973 
Arab-Israeli  wax.  In  that  war,  the  40%  of  total  Israeli 
casualties  were  caused  by  armored  vehicle  destructions.  Thus, 
tank  armor  improvement  became  the  first  priority  in  the  Israeli 
tank  corps. 

If  war  occurs  in  any  country,  the  same  problem  will  be 
faced.  That  is  the  reason  why  the  modification  of  the  conven¬ 
tional  tank  armor  is  needed  against  modern  anti-tank  weapons. 

To  improve  the  armor  protection,  the  material  properties, 
projectile  characteristics,  and  explosion  phenomena  must  be 
understood  well.  This  is  a  complicated  and  difficult  problem 
and  is  originally  a  nonlinear  problem. 

First  the  shock  wave  from  an  explosion  imposes  a  dynamic 
load  on  any  object.  That  dynamic  load  is  characterized  by 
rapidly  reached  peak  overpressure  which  then  decreases  as  the 
shock  wave  decays.  The  net  effects  of  the  load  depend  both 
on  the  nature  of  the  shock  wave  and  the  geometry  and  construc¬ 
tion  of  the  object.  For  example,  if  the  explosion  occurs  on 
thin  skirting  plate,  perforation  takes  place  by  impact  load 
and  their  impulse  loading  for  main  armor  would  be  reduced 
significantly.  The  typical  air  shock  wave  was  characterized 
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as  follows;  the  shock  wave  is  supersonic  and  generates  the 
abrupt  increasing  overpressure. 

The  non-linear  elasto-plastic  responses  of  a  clamped  flat 
plate  to  a  typical  air  shock  wave  were  investigated.  The 
nonlinear  effects  to  the  plate  responses  due  to  the  material 
and  geometric  nonlinearity  were  studied. 

In  this  study,  (1)  the  necessity  of  the  modification  of 
old  armored  vehicles  was  reviewed,  (2)  the  NASTRAN  code  was 
employed,  (3)  the  theoretical  background  of  the  nonlinear 
transient  analysis  was  described,  (4)  a  step  by  step  proce¬ 
dure  of  analyzing  the  dynamic  load  problem  by  shock  wave 
using  NASTRAN  code  was  developed,  and  (5)  sensitivity  analy¬ 
ses  were  performed  and  also  difficulties  associated  with  the 
nonlinear  analysis  were  described. 


II.  LITERATURE  REVIEW 


A.  MODERN  PROJECTILES 

An  important  version  of  the  AP  (Armor-Piercing)  round  is 
that  known  as  APDS  (Armour  Piercing  Discarding  Sabot) .  In 
this  round  the  armour-piercing  element  is  a  hard  core  of 
calibre  significantly  smaller  than  that  of  the  complete 
projectile  and  shaped  for  maximum  anti-armour  effectiveness. 
This  is  surrounded  by  a  much  lighter  casing,  or  sabot,  which 
breaks  up  as  the  projectile  leaves  the  gun  barrel.  The 
point  of  this  arrangement  is  that  the  core  can  have  a  shape 
which  differs  considerably  from  that  of  the  complete  projec¬ 
tile  and,  not  having  been  scored  by  the  rifling  in  the  barrel, 
is  aerodynaiaically  clean.  Moreover,  since  it  is  very  much 
heavier  than  the  sabot,  the  core  acquire  a  greater  energy 
than  -:ould  be  imparted  to  it  if  it  were  fired,  without  sabot, 
by  a  gun  of  appropriate  calibre  and  the  same  chamber  pressure 
as  that  of  the  larger  gun.  A  controversial  refinement 
involves  the  use  of  depleted  uranium  in  the  core.  This 
metal  has  an  exothermic  reaction  with  steel  which  aids 
penetration;  it  is  controversial  because  of  its  connection 
with  the  other  isotopes  of  uranium  and  their  radioactive 
properties . 

High  explosive  (HE)  rounds  with  contact  fuzes  can  also 
be  used  against  tanks  and  may,  if  they  are  large  enough,  be 
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effective  against  exposed  portions  of  the  tank  mechanism 
(such  as  the  tracks)  but  will  not  in  general  be  effective 
against  heavy  armour.  Two  modern  refinements  of  the  HE 
round  may  be  effectively  used  in  some  circumstances:  they 
are  the  HEAP  (High-Explosive  Armour  Piercing)  and  HESH 
(High-Explosive  Squash-Head)  rounds. 

The  HEAP  round  has  a  projectile  of  relatively  high 
velocity,  a  hardened  case  and  a  high-explosive  charge  which 
is  contact  fuzed  with  a  very  short  delay.  The  effect  of  the 
arrangement  is  that  the  projectile  either  embeds  itself  in 
the  armour  or  passes  through  it  and  then  explodes  with  an 
effect  which  is  more  damaging  than  that  of  the  simple  HE 
round. 

HESH  rounds  can  be  used  with  lower-velocity  weapons,  and 
are  so  designed  that  the  nose  of  the  shell  crumples  on  impact 
to  make  a  large  contact  area  before  the  explosive  is  detonated. 
Such  rounds  will  not,  in  general,  pierce  armour,  but  the 
explosion  sets  up  shock  waves  in  the  armour  which  in  them¬ 
selves  can  incapacitate  the  crew  and  may  also  cause  pieces  of 
metal  to  break  away  from  the  surface  remote  from  the  explosion. 
Such  fragments  fly  off  at  high  velocity  and  can  do  great 
damage  inside  a  tank  turret. 

Most  widely  used  in  infantry  weapons,  however,  is  the 
projectile  (which  may  be  a  shell,  a  rocket,  a  grenade  or  a 
guided  missile)  with  a  shaped-charge  warhead  frequently 
described  as  a  HEAT  (High  Explosive  Anti  Tank)  round.  This 


moving  at  a  speed  of  over  20,000  m/s  and  at  a  great  tempe¬ 
rature.  The  jet  takes  a  finite  distance  to  form,  hence  the 
importance  of  having  the  igniter  in  front  and  the  distance 
required  to  form  a  jet  is  known  as  the  "stand  off".  For  most 
modem  warheads  the  stand  off  distance  needs  to  be  between  2 
and  3  diameters  of  the  warhead  and  to  provide  the  full  dis¬ 
tance  sometimes  interferes  with  the  ballistics  of  the 
projectile. 

The  effects  of  a  hollow  charge  on  the  interior  of  a 
vehicle  can  be  variable.  A  very  small  jet  may  well  not 
cause  much  damage  apart  from  making  a  hole  in  the  armour  and 
destroying  whatever  is  in  its  direct  path,  but  it  is  possible 
to  "turn"  warheads  so  that  the  penetrative  effect  is  lessened 
but  there  is  more  energy  left  after  passing  through  the 
armour  and  this  energy  can  be  spread  out  into  more  of  a  fan¬ 
shaped  exit  blast  with  appropriately  enhanced  damage.  All 
hollow  charge  jets  carry  a  certain  amount  of  molten  metal 
into  the  interior  and  some  research  has  gone  into  finding 
ways  of  carrying  more  through  the  hole. 

Two  examples  of  the  Soviet  anti-tank  guided  missiles  are 
described  below.  [Ref.  1] 

1.  Sagger  Anti-Tank  Guided  Missile 

Sagger  is  the  NATO  code-name  for  the  Soviet  anti¬ 
tank  missile  which  is  also  known  by  the  US  alphanumeric 
designation  AT-3,  which  is  now  extensively  deployed  in 
Europe  and  elsewhere.  Currently  it  is  known  to  be  used  in 
at  least  three  ways: 


a)  As  a  man-portable  missile 

b)  Mounted  on  a  bracket 

c)  Mounted  on  a  vehicle 


The  characteristics  of  sagger  missile  is  as  follow. 


a) 

Type:  Wire-guided  MCLOS 

(Manual  Command  to  Line 

of  Sight) 

b) 

Weight  at  launch:  11.3 

kg 

c) 

Diameter:  120  mm 

d) 

Warhead:  HEAT 

e) 

Range:  500-3000  m 

f) 

Penetration:  400  mro 

2.  AT-4  Anti-Tank  Guided  Missile 

There  have  been  reports  for  some  years  of  a  Soviet 
SACLOS  (Semi-Automatic  Command  to  Line  of  Sight)  system,  but 
it  was  not  until  1980  that  photographs  and  a  little  informa¬ 
tion  appeared  in  some  Warsaw  Pact  military  journals. 
Characteristics : 

a)  Type:  Wire  guided  SACLOS  system 

b)  Weight  at  launch:  40.4  kg 

c)  Diameter:  120  mm 

d)  Warhead :  HEAT 

e)  Range:  2000  -  2500  m 

f)  Penetration:  500-600  mm 


Figure  2.2.  Sectioned  View  of  Sagger  Missile. 

Which  of  these  various  types  of  projectile  will  be  most 
effective  against  armour  which  depends  on  the  design  and 
thickness  of  the  armour.  Against  solid  armour  the  HEAT 
round  is  particularly  effective;  but  in  certain  circumstances 
a  sandwich  construction  is  likely  to  be  more  vulnerable  to  an 
APOS  round.  There  is  therefore  a  good  deal  to  be  said  for 
having  a  mix  of  projectiles  available  for  satisfactory 
defense  against  tanks. 

B.  CONFIGURATION  OF  ARMORED  VEHICLES 
1.  Armor 

We  can  define  the  military  armor  as  a  material 
intended  to  protect  men  and  installations  of  equipment. 
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Aluminum,  titanium  and  magnesium  alloy,  plastic  laminates, 
nylon  fabrics  and  masticstone  composites  have  all  been  used, 
but  the  most  common  form  of  armor  is  heat-treated  alloy  steel. 
This  comes  mainly  in  two  forms:  homogeneous  and  face-hardened. 

The  first  has  constant  ballistic  and  mechanical  quali¬ 
ties  throughout,  which  in  modern  tank  armor  means  high  tensile 
strength  (of  the  order  of  150,000  lb  sq  in)  and  medium  hard¬ 
ness  (i.e.  Brinnell  N  300) ;  face-hardened  armor  consists  of 
a  very  hard  first  layer  which  resists  penetration  but  lacks 
toughness  and  a  more  elastic  core  to  support  the  face.  Naval 
deck  armor  and  modern  tank  armor  are  homogeneous,  while  face- 
hardened  armor  is  still  used  for  some  light  armored  vehicles 
and  side  protection  for  naval  vessels. 

Steel  armour  usually  means  rolled  or  wrought  plates, 
but  castings  are  used  for  tank  turrets  and  even  for  entire 
tank  hulls,  as  in  the  M.48  and  M.60  tanks  (Steel  armor  is 
usually  a  nickel-chrome-molybdenum  alloy,  but  other  alloys 
including  roanganese-silicon-chromium-nickel  steel  have  been 
used,  especially  in  Soviet  tank  armor.). 

Typically,  a  well-designed  hollow-charge  warhead  can 
be  expected  to  penetrate  armor  up  to  five  times  the  diameter 
of  the  warhead.  The  best  defense  against  such  attack  is  not 
to  make  thicker  armor,  but  to  dissipate  the  jet  before  it 
reaches  the  surface.  To  do  this  a  light  screen  is  carried 
some  distance  in  front  of  the  main  armor,  of  sufficient 
thickness  to  set  off  the  fuze  of  any  hollow-charge  projectile. 
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Quite  naturally  such  screens  do  not  last  long,  but  they  are 
easily  replaced  and  provide  almost  complete  protection  for 
the  armor  behind.  A  common  use  of  such  a  technique  is  the 
thin  skirting  plates  carried  over  the  tracks  of  all  modern 
tanks.  Before  leaving  this  subject  it  is  appropriate  to 
mention  a  recent  British  armor  development,  the  •Chobham' 


Figure  2.3.  The  Effect  of  Anti-Tank  Projectiles  on  "Cobham" 
Armor. 

armor.  This  is  claimed  to  offer  greatly  enhanced  resistance 
to  all  forms  of  projectile  as  is  indicated  by  the  accompanying 
illustrations.  Details  of  the  construction  of  this  armor  have 
not  been  released,  but  it  seems  probable  that  it  is  a  sand¬ 
wich  made  up  of  two  sheets  of  armor  plate  separated  by  a 
highly  dispersive,  possibly  granular,  filling.  The  disad¬ 
vantage  of  this  armor  lies  in  its  bulk.  To  gain  the  protec¬ 
tion,  it  has  to  be  thick  and  fairly  heavy;  tanks  are  already 
large  enough  and  the  addition  of  armor  plating,  which  is  200 
mm  or  more  in  thickness  makes  a  startling  difference  to  the 
outline  and  bulk  of  the  entire  vehicle. 


2.  Current  Modification  of  Old  Tanks 

Determining  which  one  should  we  choose,  a  new  one 
or  modification  of  an  old  one  usually  depends  on  cost  and 
political  situations. 

The  U.S.  completed  conversion  of  M48  tanks  to  M48A5 
standards  in  1980.  This  conversion  basically  included  adding 
diesel  engines  to  M48Als  and  adding  105-mm  guns  to  the  tur¬ 
rets  of  M48Als  and  M48A3s.  Numerous  minor  changes  were  also 
made.  But  little  improvement  in  armor  protection  has  been 
made.  If  we  want  to  take  more  modifications  for  M48A3s  and 
M48A5s  tanks,  a  spaced  applique  armor  should  be  added  and 
include  a  1,200  hp  engine,  improved  transmission,  better 
fire  control,  and  a  hydropneumatic  suspension  system. 

The  following  example  shows  modern  modifications. 

The  High  Performance  M60  MBT  has  been  developed  as  a 
private  venture  by  the  General  Products  Division  of  Teledyne 
Continental  Motors.  It  is  essentially  an  M60  series  MBT 
fitted  with  a  new  powerpack,  new  suspension  and,  as  an  option, 
additional  armor.  These  modifications  can  be  carried  out  on 
existing  M60  series.  Teledyne  believes  that  to  fit  an  M60A1 
with  the  new  powerpack  would  cost  under  $200,000,  while  the 
added  armor  would  cost  just  $40,000. 


Israeli  experience  in  the  1967  campaign  proved  that 
mobility  was  not  substitute  for  armor  protection  and  they 
therefore  decided  at  an  early  stage  that  the  main  emphasis 
would  be  placed  on  armor,  with  firepower  and  mobility  second 
and  third  priorities.  After  that  they  spent  a  lot  of  money 
to  figure  out  this  protection  problem.  Total  cost  of 
research,  development,  trials  and  the  building  of  the  proto¬ 
type  vehicles  was  about  $65  million. 

In  May  1977  Israel  finally  announced  that  it  had 
developed  a  new  MBT  called  the  Merkava  which  was  made  of 
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cast  and  welded  armor  with  a  well-shaped  glacis  plate. 

Behind  the  first  layer  of  cast  armor  is  a  space  filled 
with  diesel  fuel  and  then  another  layer  of  armor.  This 
spaced  armor  gives  the  tank  protection  from  HEAT  projec¬ 
tiles  and  ATGWs  (anti-tank  guided  weapons) . 

They  also  modified  their  old  tanks  such  as  M60s, 
Centurions ,  etc . 

The  add-on  and  readily  detachable  protective  arrays 
seen  on  Israeli  M60s  and  Centurions  contain  ceramic  tiles, 
would  also  mark  a  first  in  terms  of  operational  deployment. 

The  possibility  cannot  be  excluded  that  armor  of  both  kinds 
is  involved. 

While  no  answers  are  to  be  found  in  external  appear¬ 
ance,  the  configurations,  as  seen  in  the  accompanying  Defence 
Attache  photographs,  show  obvious  differences  in  the  thick¬ 
ness  and  shape  of  the  packs  provided  for  M60  and  Centurion 
respectively.  The  Centurion's  additional  armor  is  mainly  in 
the  form  of  quite  shallow  panels  fitted  to  the  brow  of  the 
turret  and  thicker  pieces  to  the  glacis  and  the  mantlet 
either  side  of  the  main  gun;  a  hole  is  included  in  the  left- 
hand  piece  to  permit  the  firing  of  the  coaxial  machine  gun. 
Wedge-shaped  packs  are  mounted  on  top  of  the  storage  boxes 
fitted  forward  above  the  track  guards. 


the  turret.  Here  the  packs  are  thick  and  vertically  deep, 
with  smaller  and  thinner  sections  patched  in  —  with  no 
great  concern  for  aesthetics  —  to  the  upper  surfaces  which, 
because  of  the  turret's  curvature  and  varying  angles  to  the 
perpendicular,  not  to  mention  obstructions  such  as  lifting 
eyes,  cannot  be  covered  by  the  main  boxes.  As  with  the 
Centurion,  panels  are  fitted  to  the  front  of  the  turret  roof 
and  to  the  glacis;  in  the  latter  case,  because  the  M60's 
basic  glacis  is  curved,  it  has  been  necessary  to  fit  a 
straightening  frame  to  which  the  packs  can  be  attached.  The 
M60's  installation  is  completed  by  the  fitting  of  wedge- 
shaped  or  angled  sections  either  side  of  the  driver's  hatch 
to  help  protect  the  vulnerable  area  where  turret  meets  hull, 
with  further  angled  sections  each  side  of,  and  thinner  plates 
above,  the  mantlet. 

Although  it  can  be  stated  with  a  fair  degree  of 
certainty  that  the  add-on  arrays  are  intended  to  help  defeat 
shaped-charge  high-explosive  anti-tank  (HEAT)  attack,  any 
discussion  of  what  the  packs  and  panels  consist  of  must,  in 
the  absence  of  any  word  of  guidance  from  the  Israelis  -  who 
on  the  contrary  are  happily  watching,  if  not  promoting,  the 
circulation  of  a  number  of  often  fanciful  and  frequently 
mutually  exclusive  'solutions'  to  the  puzzle-fall  short  of 
a  definitive  answer. 

While  one  source  suggested  that  the  boxes  were  empty 
(and  the  large  ones  on  the  M60  do  sound  quite  hollow  when 


thumped)  the  probability  is  that  there  is  more  to  it  than 
that.  In  view  of  solutions  that  have  been  adopted  else¬ 
where,  though  usually  in  the  context  of  new-build  tanks,  it 
seems  quite  likely  that  the  Israelis  have  developed  their  own 
ceramic  tiles  and  that  the  boxes  are  partially  filled  with 
them.  The  alternative  of  a  form  of  active  armor  is  suggested 
more  by  the  vehemence  with  which  several  sources  in  Israel 
stated  this  to  be  the  case  than  by  any  objective  evidence. 

A  number  of  different  types  and  configurations  of 
active  armor  are  possible.  An  example  is  a  sheet  of  explosive- 
laid  of  top  of  a  steel  plate,  with  a  covering  for  environmental 
protection  on  top  -  which  detonates  when  hit  by  the  incipient 
jet  of  a  shaped-charge  warhead  and  in  the  process  disrupts  the 
jet's  further  formation.  If  the  Israelis  are  indeed  deploying 
a  form  of  active  armor,  such  a  fact  would  carry  the  implica¬ 
tion  that  they  have  sufficiently  overcome  the  formidable 
difficulties  that  seem  to  have  prevented  larger  and  better- 
endowed  countries  from  doing  likewise;  the  United  States  and 
Soviet  Union  are  among  the  nations  that  are  known  to  have  been 
engaged  in  related  research  and  experimentation  for  many  years. 

C.  ELEMENTARY  DETONATION  THEORY 

1.  The  following  assumption  will  be  taken  for  the 
simplest  analysis  of  a  detonation:  [Ref.  2  ] 


a)  The  flow  (motion  of  the  shock  front)  is  one 


dimensional 

b)  The  plane  detonation  front  is  a  jump  discontinuity 
at  which  chemical  reaction  is  assumed  to  be 
instantaneous . 

c)  The  discontinuity  is  time-independent,  so  that  the 
product  composition  (at  least  directly  behind  the 
shock  front)  is  likewise  time-independent. 

u  represents  the  particle  velocity  of  blast  wind,  which  is 
that  of  the  products  behind  the  shock  front.  Subscript  x 
represents  the  initial  values  of  the  properties,  y  those 
after  shock.  D  is  the  detonation  speed  and  p  is  the  shock 
wave  overpressure. 

From  mass  conservation  across  the  front  we  can  write 


P  D  =  P„  (D-u  ) 


Py*py 


Figure  2.7.  Explosive  Shock  Front 
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and  from .conservation  of  momentum, 

Py  -  Px  -  Px  upD  (2) 

Eliminating  up  from  the  two  equations  we  obtain  the  equation 
for  the  Rayleigh  line: 

PX2I>2  -  (py  -  px)/(vx  -  vy)  -  0  (3) 

(Observe  that  in  part  of  the  previous  equation  we  replace  1/p 
by  the  specific  volume,  v) .  Since  for  a  detonation,  py  >  px, 
it  follows  that  vy  <  vx  (for  real  values  of  the  detonation 
velocity) ;  and,  indeed,  it  is  observed  as  a  rule  that  there 
is  about  a  one-third  increase  in  density  on  detonation.  If 
this  condition  is  substituted  in  eq  (1) ,  it  follows  that 
0  <  ^  <  d,  so  that  the  products  travel  toward  the  detonation 
front,  though  trailing  behind. 

If  the  Rayleigh  equation  is  solved  for  py  we  obtain 
Py  -  <PX  +  PXD2)  -  Px2D2  vy  (4) 

wherein  py  is  linear  with  vy,  the  slope  of  the  line  being 
2  2 

-p  D  .  The  conservation  conditions  thus  require  that  the 
«■» 

initial  and  final  states  of  a  detonation  must  lie  at  the 

termini  of  the  Rayleigh  line. 

Conservation  of  energy  leads  to 

E ( P  , v  ,  A=l) +P  v  +1/2 (D-U  ) 2  =  £ (P  ,1/x, A=0) +P  v  +1/2D2  (5) 

y  y  yy  p  *  A  A 

in  which  A  =  the  extent  of  chemical  reaction  (A  =0  for  no 
reaction  and  1  for  complete  reaction) .  Eliminating  u  and 
D  from  Eq  (5),  suing  Eqs  (1)  and  (2)  leads  to  the  Hugoniot 
curve : 
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E(WA=1)  "  (PX'VX'A=0)  “  1/2  tPy+V  {Vx'Vy)=0  (6) 

A  EH  =  1/2  (Py+Px)  Cvx-V  )  (7) 

AEH  is  the  internal  energy  of  the  detonation  products  at  the 
very  high  pressure  and  temperature  of  the  detonation  state 
less  the  internal  energy  of  the  unreacted  explosive  at  ambient 
temperature  and  pressure.  It  may  be  written  as  the  sum  of  two 


terms : 


AEh  =  AEt  +  AEp  (8) 

where  AE^,  is  the  internal  energy  change  for  the  conversion 
of  reactants  to  products  isothermally ,  and  AE  is  that 

Br 

required  to  warm  and  compress  the  products  from  ambient  p 
and  T  to  the  detonation  values.  AET  is  often  symbolized 
as  -Q  in  explosive  literature  (Q  being  the  heat  evolved  in  an 
explosion  when  corrected  to  isothermal  conditions) .  If  the 
products  were  ideal  gases ,  then  AE  would  be  given  |c  dT, 

r  " 

and  this  approximation  is  often  made,  though  the  products  at 
pressures  of  hundreds  of  bars,  are  far  from  being  ideal  gases. 

2.  The  Hugoniot  curve  is  that  of  a  hyperbola  in  the  p-v 
plane.  For  the  non-reacting  case  (\=0)  the  hyperbola  passes 
through  the  initial  state,  p  ,v  ;  but  for  the  case  we  are 

A  A 

considering  here  (X=l)  it  does  not.  The  conservation  laws 
require  that  the  final  state,  however,  lies  on  the  Hugoniot 
curve  as  well  as  on  the  Rayleigh  line.  We  consider  below 


several  cases: 
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For  a  too  small  D^c^)  the  Rayleigh  line  does  not 
intersect  the  Hugoniot  and  thus  no  solution  exists.  For 


Figure  2.8.  The  Hugoniot  Curve. 

large  D  (=0^) ,  there  are  two  intersections,  called  strong 

(labelled  S)  and  weak  (W) .  For  one  particular  value  of 

D  (=D  .)  the  Rayleigh  line  is  tangent  to  the  Hugoniot. 

CJ 

It  can  be  shown  that  the  relative  particle  velocity, 

D-u  ,  is  subsonic  at  S  (D=u  <  a  ) ,  supersonic  at  W  and 
p  p  . 

sonic  at  the  point  of  tangency,  called  the  Chapman-Jouguet 

point:  (D-u  )  .  »  a  (sonic  velocity  in  the  shocked  medium) . 

P  c)  y 

There  exist  a  number  of  computer  codes  based  on  such 
non-ideal  gas  treatments  (RUBY  and  TIGER  are  two  of  these)  and 
moderate  success  has  been  achieved  in  computing  detonation 
properties  from  these  codes.  The  equations  are: 


(9) 


P  =  kp  2  4> 
y  Hx  y 

D  =  A*3*  (1  +  BPx)  (10) 

<t>  =  nmV5  (11) 

with  k  =  0.762,  A  =  22.3,  B  =  0.0013  (SI  units).  Q,  the 
heat  of  explosion  in  J/kg;  N,  the  number  of  moles  of  gas  per 
kg;  and  M,  the  average  molar  mass  of  the  gas  (kg/mole  are  all 
found  from  an  assumed  stoichiometry  of  decomposition,  combined 
with  thermochemical  principles.  Although  N ,  M,  and  Q  are  all 
strongly  affected  by  the  assumption  made  as  to  stoichiometry, 
it  turns  out  that  the  important  parameter,  <p,  is  rather 
insensitive  to  the  assumption. 


III.  ANALYSIS  OF  DYNAMIC  RESPONSE  OF  PLATE  BY  SHOCK  WAVES 


A.  THEORETICAL  BACKGROUND 


1.  Nonlinear  Static  Analysis 


The  fundamental  method  used  to  obtain  a  static  solu¬ 


tion  is  to  minimize  the  error  vector  {S},  given  by  the 
following  equation:  [Ref.  7] 


{&}  -  {p}+  (q}-{f } 


(12) 


where  {6}  is  the  error  vector  of  unbalanced  forces  acting 
at  all  grid  points,  {P}  is  the  vector  of  applied  external 
loads,  which  may  change  with  displacements,  (Q }  is  the  unknown 
vector  of  constraint  forces  due  to  single  and  multipoint  con¬ 
straints,  and  {F}  is  the  vector  of  grid  point  forces  generated 
by  element  motion  and  stress.  The  terms  are  functions  of 
displacement,  temperature,  and  stress  history. 

It  must  be  noted  that  degrees  of  freedom  not  involved 
in  constraints  produce  null  terms  in  {Q}  and  that  dependent 
constraint  points  produce  no  errors.  Therefore,  when  Eq  (12) 
is  reduced  to  the  solution  coordinates  the  constraint  forces, 
{Q},  disappear.  Terms  in  the  vector  {F}  are  dependent  on  the 
deformations  of  the  finite  elements.  In  linear  static  analysis 
the  vector  {F}  becomes 

{Flinear}  *  (kPHu>  *  <PT<T>  >  U3) 
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where  [k^]  =  The  linear  stiffness  matrix 
T 

{P  }  =  The  "thermal  load"  vector,  i.e.,  a 

vector  of  grid  point  which  produces  the 
same  element  strains  as  the  temperature 
distribution, 

{u}  =  The  vector  of  grid  point  displacements. 

To  obtain  a  nonlinear  solution,  a  "tangent"  matrix 
[K]  is  used  in  which  the  terms  are  derivatives  in  the  form 

6Fi 

Kij  *  ^<$U .  ^  u  =  u 
J  :  r 

.The  matrix  [K]  contains  both  geometric  and  material  non¬ 
linear  effects.  The  approximation  occurs  because  the  forces 
are  highly  variable  with  displacements  and  nonsymmetric  terms 
are  approximated.  If  {u1}  is  a  known  displacement  which 

produces  an  error  {5(ux)},  then  Newton's  method  may  be  used 

i+1 

to  predict  a  new  vector  {u  }  with  a  smaller  error.  For 
a  displacement  {u}  near  a  known  solution,  {ux>,  the  non¬ 
linear  force  is  approximately 

{F (u) }  %  {F (ux) }  +  fK] {u-u1}  (15) 

The  substitution  of  Eqn  (15)  into  Eqn  (12)  with  {5}  =  0 
yields  the  result 

[KHu-U1}  ^  {P}  +  {Q}  -  {F  (uX)  }  (16) 

which  can  be  written  in  terms  of  the  error  {5}  to  provide 
the  format  of  the  Newton-Raphson  iteration  method: 

[KHu^-u1}  =  Ufu1)}  (17) 
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When  [K]  may  be  inverted,  Eqn  (17)  may  be  used  to  solve 
for  {ui+1}  =  {(u)},  a  new  estimate  of  the  solution. 

Equation  (17)  provides  an  incremental  form  for  solution 
iteration.  Another  form  of  the  iteration  equation  is 
obtained  from  Eqns  (12)  and  (17) : 

[K] {u1+1-u1} = {pQ  +  pT}  +  {Q}  -  {f1}  -  [K]  {u1}  (18) 

where 

{f1}  =  {Ffu1)}  -  [K]  {u1}  +  {pT}  -  Ap (u1)  (19) 

[K1]  is  the  nonlinear  "reference"  matrix,  available  from  the 
elements  and  Aplu1)  is  the  change  in  loads  due  to  grid 
point  motions.  The  advantage  to  this  form  is  that  the  second 
order  "corrective"  load  vector  {f1}  equals  zero  when  the 
forces  {F>  are  linear  and  when  [K]  contains  linear  elastic 
terms  (see  Eqn  13) .  Because  the  { f }  forces  are  calculated 
at  the  element  level,  all  known  linear  elements  and  super¬ 
elements  may  be  bypassed  in  their  calculation.  This  is  very 
efficient  when  a  major  portion  of  the  structure  remains  linear 
However,  the  advantages  for  using  Eqn  (17)  are  that  the  load 
error  and  incremental  displacements  (u^+^  -  u^}  are  directly 
available  for  testing  convergence.  Therefore,  in  MSC/NASTRAN , 
both  Eqns  (16)  and  (17)  are  used  internally. 

In  order  to  avoid  calculating  the  vector  {FCu1)}, 

Eqn  (19)  is  solved  for  {F(ux)}  and  substituted  into  Eqn  (12) 
to  obtain 

{di}  =  {5( u1)}  =  {p  +  pT}  -  {f1}  -  [Kr]  (u1)  (20) 
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Subtracting  {61}  for  two  successive  values  of  i  in 
Egn  (20)  to  eliminate  the  constants  a  simple  equation  is  ob¬ 
tained  for  the  load  error  in  the  form 

{51}  =  {S1"1}  -{f1  -  fi_1}  -  [k]  (u1  -  u1'1}  (21) 

This  form  provides  a  simple  iteration  procedure,  namely 

a)  For  the  first  iteration  after  a  change  in  the  load  vector, 
use  Eqn  (17)  and  (20) .  Note  that  the  initial  values  for 
starting  from  rest  are  {f°}  =  {u°>  =  0:  calculate 

{6°}  -  (p>  +  {pT>  and  solve  [K]  {u1}  =  {<5°}  (22)  (23) 

b)  then  compute  f^  directly  from  the  element  and  load 
tables: 

{f1}  =  FCu1}  -  [Kr]  {u1}  +  {pT}  -  (Ap (u1) }  (24) 

c)  then  compute  {51}  from  Eqn  (21) . 

{51}  *  5°  -  {f1}  +  {f°}  -  [Kr]  {u1  -  u°}  (25) 

where  {f°}  =  0  by  assumption. 

2 

d)  Use  Eqn  (17)  with  i=l  again  to  compute  (u  } 

[K]{u2  -  u1}  •  SCu1}  (26) 

e)  Repeat  steps  b,  c  and  d  to  convergence,  incrementing  the 
superscripts  by  one  at  each  iteration  loop. 

2.  Nonlinear  Transient  Analysis  Using  MSC/NASTRAN 

SOL  99  which  is  a  rigid  format  of  MSC/NASTRAN  code 
for  nonlinear  transient  problem,  is  the  dynamic  comparison  to 
nonlinear  static  systems,  providing  for  combined  material  and 
geometric  nonlinear  analysis  in  the  time  domain.  SOL  66  which 
is  a  rigid  format  for  nonlinear  static  problem  and  SOL  99 
share  the  same  element  computer  code  and  solution  techniques 
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and  provide  similar  data  storage  and  restart  facilities. 

The  transient  solution  is  performed  in  a  stepwise  manner  with 
time  steps  replacing  load  steps  and  with  effects  of  mass, 
damping  and  unsymmetric  matrices  added  to  the  stiffness 
matrices . 

Because  the  calculations  are  identical,  many  element 
and  material  discussions  in  the  previous  sections  also  apply 
to  the  transient  analysis  system.  The  finite  elements  auto¬ 
matically  provide  linear  inertial  forces  and  structural 
damping,  based  upon  the  initial  state,  as  well  as  nonlinear 
forces  due  to  geometry  and  material  effects. 

We  will  develop  the  basic  equations  for  the  transient 
solution.  The  methods  were  chosen  to  be  compatible  with 
existing  MSC/NASTRAN  static  nonlinear  solutions  as  well  as 
with  the  existing  linear  transient  methods.  Thus,  the 
Newmark  Beta  method  for  transient  integration  is  combined 
with  the  modifications  to  Newton's  method  for  nonlinear 
solutions.  The  additional  iteration  steps  provide  equilib¬ 
rium  solutions  at  each  time  step,  thereby  guaranteeing 
stability  and  accuracy  for  arbitrary  time  step  size.  The 
basic  equations  are  developed  below,  followed  by  a  discus¬ 
sion  of  the  stability  limits. 

a.  Basic  Equation 

Adding  dynamic  loads  to  the  basic  equation  given 
previously,  we  obtain  a  load  equilibrium  error  vector  {6n> 
at  time  step  n  by  equation: 


where  {Pn>=  "Avera9eH  l°ad  over  the  time  period 
<fcn-l  "  fcn  "  W 

•  •  • 

{Un}  {Un)  =  Corresponding  acceleration  and  velocity  vectors. 

{Fn>=  "Average"  elasto-plastic  element  total  force 
vector. 

The  above  equation  is  solved  at  the  "reduced" 

(u^)  displacement  vector  size.  The  approximation  errors  due 
to  dynamic  reduction  methods  are  not  included  in  the  error 
vector  {<5n}.  Applying  Newmark's  averaging  method  over  a 
finite  time  period,  fcn=i  <  t  <  tR+^  the  "static"  forces  are 

{Fn>  «  { |3F(un+^)  +  (l-26)F(un)  +  gF(un-1)}  (28) 

where  $  *  Newmark  Beta  operator  and  F(un)  is  the  nonlinear 
force  due  to  a  generalized  displacement  vector  {Un> .  An 
identical  definition  occurs  for  {pn>  from  the  applied  loads 
at  each  time  step. 

From  central  finite  differences,  the  acceleration 
and  velocity  vectors  are: 

{un}  =  772  {Aun+l  ”  *un}  =  77?  {dn}  {29) 

At  At 


where  {Au. 


tun>  *  Tit  Uun+1  +  AV 

{un+l  -  un> 


for  small  angles  and  At  =  fcn+i  ”  fcn  is  the  "time  step  size". 
For  large  angle  changes  we  require  that  (u),  {u},  and  {Au} 
are  vectors  in  the  global  coordinate  direction.  The  vector  {u} 


contains  gimbal  angles  requiring  a  tranformation  R(un) 
such  that 


{lW  =  lun>  +  tRn}  {ilW 


(31) 


the  load  error  in  terms  of  the  current  vectors  estimates: 

'V  =  ‘V  -  M  +  2St  B'  {dn>  -1st-  B1UV 


(S£n+1  +  (1-28)  Fn  +  BF^} 


(32) 


and 


lAun+l  }  -  {4un  +  dn} 


(33) 


At  any  time  step  the  vectors  {P  },  (AU  },  {F  }  ,  and  (F  .} 

n  n  n  n-i 

are  known.  In  general  the  vectors  {dn>,  (un+1>  and  fn+l  = 
F(un+^)  must  be  found  either  by  approximating  {Fn+^}  or  by 
an  iterative  search. 

Note  that  for  a  linear  solution  (Fn>  equals 
[K]{u  },  where  [K]  is  the  reduced  stiffness  matrix.  In  this 
case  Eqns  (32)  and  (33)  may  be  used  to  solve  for  {un+1> 
directly  by  setting  {<$n}  =  0.  With  8=1/3  this  will  produce 
the  linear  SOL  27  MSC/NASTRAN  results. 

Nonlinear  Iterations 

For  a  nonlinear  solution,  Equns  (32)  and  (33)  may  be 
solved  with  Newton's  method  (or  a  modified  version) .  Using 
{dn>  as  the  primary  solution  variable  provides  the  following 
iteration  algorithm. 


First  assume  a  linear  approximation: 


{$  1+1)  ~  {5  }  ttt-  (d 1+1  -  d  L}  =  0 
n  n  da  n  n 

n 


{A,}  -  35n 


Then,  the  estimated  displacement  change  is 

{dn+1  -  dn>  -  (Ai^t  Vdn>  1  <3«> 

The  algorithm  is  identical  to  the  static  case  except  that 
instead  of  a  tangent  stiffness  matrix,  the  "left  hand  side" 
matrix,  obtained  from  Eqns  (32)  and  (34)  is 


{A  }=  [Ay  M  +  A-  B  +  6K} 


3F  , 


where  [K]  ^  is  the  current  tangent  stiffness  matrix. 


At  a  new  time  step  the  load  iterations  may  be  started  by 

O 

assuming  that  {AUn+^}  =  (AUn)  anc*  therefore: 


{Sn}  ={0} 


However,  to  be  consistent  with  displacements,  we  must  extra¬ 
polate  the  nonlinear  forces: 

'Cl1  S  (Pn  +  <Fn  *  ^-l11  '  (2Pn  *  F„-l>  <3 

Substituting  Eqns  (38)  and  (39)  into  Eqn  (32)  for  i=0  results 
in  the  first  estimate: 


(5°>=  {P  }  +  [A,]  { Au  }  -{F  } 
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where 


tA2J  *  t-  Bl 


(41) 


Here  A  is  a  precalculated  matrix.  {AUn>  and  (Fn>  are 
previous  results.  For  subsequent  iterations,  (i>0)  Eqn  (32) 
becomes : 


(sb  =  - 

n  n 


»<Cr 


Fn+X> 


(42) 


where  .  . 

[A3]  =  [-i*  M  +  B]  (43) 

At 

At  each  time  step  the  code  will  iterate  on  Eqns 
(36)  and  (42)  until  {S^}  passes  the  convergence  tests  or 
the  number  o£  passes  reaches  an  iteration  limit.  With  a 
single  step,  i=l,  calculating  only  6^,  the  results  will  be 
identical  to  the  existing  NOLIN  results  in  MSC/NASTRAN.  For 
faster  convergence  the  iterations  may  continue,  the  matrices 
may  be  updated,  and/or  the  ^ime  step  size  may  be  reduced. 

3 .  Factors  Determining  Response  to  Transient  Loading 
The  response  of  a  metal  to  dynamic  loading  depends 
upon  three  elements  which  are  characteristic  of  the  material: 
the  dynamic  constitutive  relation,  the  pressure-temperature 
phase  diagram,  and  the  dynamic  fracture  criteria.  In  the 
absence  of  phase  changes  or  fracture,  the  material  response 
is  governed  completely  by  the  constitutive  relation,  or 
dynamic  equation  of  state.  The  constitutive  relation  is 
characterized  as  "elastic",  "viscoelastic",  "visco-elastic- 


plastic",  etc.  It  depends  on  macroscopic  parameters  such  as 


Young's  modulus  and  the  viscosity  coefficient  which,  on  a 
microscopic  scale,  arise  from  atomic  interactions  and  dislo¬ 
cation  processes.  The  phase  diagram  determines  whether  a 
given  loading  history  will  bring  the  metal  into  a  pressure- 
temperature  regime  where  phase  transitions  are  possible.  If 
this  occurs  during  shock  loading,  multiple  shock  waves  may 
result.  Microcracks  and/or  voids  begin  to  appear  in  a  metal 
under  tension  whenever  the  criteria  for  nucleation  and  growth 
of  damage  are  met.  These  criteria  are  functions  of  material 
parameters  and  the  stress  history. 

In  this  chapter  the  dynamic  constitutive  relations 
are  discussed.  It  is  the  most  important  element  to  determine 
the  response  of  a  metal  to  dynamic  loadings. 

In  macroscopic  or  rheological  language,  iron  can  be 
described  as  a  viscoelastic-plastic  material,  i.e.,  it 
responds  elastically  until  yielding  (with  upper  and  lower 
yield  points) ,  and  is  strongly  rate-dependent. 

The  constitutive  relation,  or  dynamic  equation  of 

state,  thus  gives  the  stress  as  a  function  of  strain  and 

time,  or,  alternatively,  can  be  written  as  a  differential 

•  • 

equation  in  o,e,a,e  and  an  important  goal  of  present  day 
research  is  to  develop  the  capability  to  derive  this  macro¬ 
scopic  constitutive  relation  by  averaging  the  effects  of 
microscopic  processes . 

This  is  analogous  to  the  way  in  which  the  thermo¬ 
dynamic  equation  of  state  of  a  gas  is  derived  by  averaging 
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the  effects  of  the  molecular  collisions.  Once  developed, 

this  capability  would  be  an  extremely  powerful  tool,  since 

the  response  of  a  metal  to  dynamic  loads  could  be  predicted 

from  known  microscopic  material  properties.  Although  this 

goal  has  not  been  attained  as  yet,  there  has  been  significant 

progress  during  the  last  few  years  in  this  direction  for  the 

case  of  crystalline  materials,  in  which  dislocation  processes 

play  the  role  of  the  "molecular  collisions".  In  iron,  the 

8  2 

density  of  mobile  dislocations  is  of  the  order  of  10  /cm  , 
so  the  number  of  dislocations  present  is  certainly  large 
enough  to  make  a  statistical  approach  valid. 

The  basic  approach  is  based  on  the  relation  from 
dislocation  theory:  [Ref.  3] 

yP  =  N  bv  (44) 

m 

where  y  is  the  plastic  shear  strain  rate 
1  *?i 

(y-  •  *  s-  (3 —  +  3 — L  )  ,  where  ?.  is  the  displacement  in 

13  £  dXj  dX^  1 

the  x^  direction) ,  Nm  is  the  mobile  disclocation  density, 
b  is  the  Burgers  vector,  and  v  is  the  average  disclocation 
velocity.  This  equation  provides  the  link  between  microscopic 
dislocation  motion  and  macroscopic  plastic  straint  rate,  and 
is  based  on  the  fact  that  it  has  been  found  possible  for  many 

p 

crystalline  solids  to  express  and  v  as  functions  of  y 

and  applied  macrosopic  shear  stress,  t. 
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The  constitutive  relation  is  obtained  by  combining 
Eqn  (44)  with  the  incremental  Hooke's  law 

doij  -  +  xd4  <«> 

under  the  assumption  that  the  total  strain  is  the  sum  of 
elastic  and  plastic  components: 

£ij  =  £ij  +  £ij  (46) 

where  A  and  p  are  the  Lame  constants,  is  the 

Kronecker  delta,  and  a  repeated  index  indicates  summation. 
The  plastic  strain  is  assumed  to  cause  no  volume  change,  i.e, 

4  - 0  <47> 

The  combination  of  Eqns  (46) -(48)  yields,  after  taking  the 
time  derivative: 

“ij  "2  »*ij  -  UKK5ij  *  2  (48> 

In  order  to  combine  Eqn  (48)  with  Eqn  (44) ,  we  must 
first  specify  the  loading  geometry.  In  plate-slap  tests 
under  conditions  of  uniaxial  strain,  all  total  strain  com¬ 
ponents  except  e  in  the  direction  of  wave  propagation  are 
zero,  and  evv-  e  -  Then  equation  (48)  becomes 

i\J\  X 


ox  -  (A+2u)ex  =  -2  U  £x 


Also,  in  this  case,  the  plane  of  maximum  shear  stress  in 
isotropic  media  is  inclined  at  45°  to  the  direction  of  shock 


propagation ,  and 


3  ;p 
T  £x 
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The  dynamic  constitutive  relation  for  this  loading  geometry 
is  thus  obtained  by  combining  Eqn  (44) ,  (49) ,  and  (50) : 

ax-  (X  +  2u)ex  =  -  |  nNmbv  (51) 

If,  X,  U/  and  the  right-hand  side  of  Eqn  (51)  are  known 
as  functions  of  macroscopic  stress  and  strain,  then  the 
equation  can,  in  principle,  be  solved  for  any  problem 
involving  uniaxial  strain. 

4 .  Dynamic  Plastic  Deformation  of  Plates 

We  consider  cases  of  the  impulsive  loading  of  thin 
metal  plates,  a  simple  situation  in  which  an  exponential 
shock  wave  impinges  on  a  clamped  square  plate.  In  this  case, 
the  plastic  deformation  is  supposed  to  be  brought  about  by 
bending  only  and  combined  bending  and  tension  is  not  analy¬ 
tically  considered.  The  pressure  in  the  shock  wave  may  be 
assumed  to  be  given  by  (Ref.  4], 

P  =  Pm  exp ( -t/ x ) 

where  p  is  the  peak  pressure  at  a  point  at  the  head  of  the 

wave,  t  denotes  time  and  i  is  a  time  constant  the  time 

taken  for  the  pressure  to  fall  to  Pm/e,  which  for  high 
explosives  is  ~30  ysecs.  We  obtain  simply,  some  of  the 
results  of  Cox  and  Morland  (Ref.  3  ],  who  theoretically 
investigated  the  dynamic  deformation  of  plates,  side  length 
2a,  position-fixed  at  their  periphery,  when  subject  to  a 
uniform  constant  transverse  pressure  p  for  time  t.  Their 
detailed  examination  showed  that  two  categories  of  load  may 
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be  distinguished;  medium  load  when  p  <  p  <  2p  and  high 

s  s 

load  when  2p  <  p,  where  p  denotes  the  pressure  required 
s  s 

just  to  enforce  plastic  collapse  of  the  plate,  statically 
For  p  we  have,  see  Fig.  3-1. 

4ay'T  *  W/?  •  Mp  =  ^  •  4a2  *aw  •  Ps  (52) 


i  .e. 


6M^/a" 


For  simplicity  we  deal  only  with  their  first  case, 
and  for  the  second,  more  complicated  one,  the  reader  is 
referred  to  the  original  source.  For  this  medium  range  of 
loading  two  sub-divisions  are  required:  (1)  when  0<t<x 
and  (2)  when  t>T.  Use  is  made  of  the  results  obtained 
above  for  static  situations  using  hinge  lines,  and  we  take 
over  whatever  results  may  be  useful.  As  previously,  the 
analysis  pays  no  special  attention  to  material  rate  effects. 
Further,  it  is  implicit  in  this  approach  that  the  static 
and  dynamic  modes  of  deformation  are  identical. 


FRONT 


\  NET  ON  THE 
!l  FRAME 


Figure  3.1.  Dynamic  Loading  of  Square  Plate 
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a. 


Medium  Load:  0  <  t  <  t 


Assume  the  mode  of  deformation  of  the  plate  is 
that  of  Fig.  3.2  and  that  the  plate  diagonals  become  hinge 
lines.  Since  p  is  constant  and  greater  than  pg,  and 
because  the  four  equal  parts  of  the  plate  act  as  rigid  bodies 
each  rotating  about  its  boundary  or  side,  they  must  undergo 
an  angular  acceleration,  a.  If  changes  in  plate  geometry 
throughout  the  deformation  are  neglected,  the  resistance  to 
deformation,  which  arises  from  the  yielding  taking  place  in 
the  hinge  lines,  will  be  constant;  this  resistance  is  related 
to  p  ,  the  pressure  to  cause  static  yield.  Thus  the  excess 
load  is  constant  and  hence  a  must  be  constant.  The  work 
done  by  the  excess  load  will  manifest  itself  as  kinetic 
energy  acquired  by  the  four  rigied  regions  of  the  plate,  i.e. 
4  *  j  Iui2,  where  I  is  the  moment  of  inertia  of  a  triangle 
about  a  side  of  the  plate  and  u  is  its  current  angular 


velocity.  Because  a  is  constant,  at  the  end  of  time  the 
plate  centre  will  have  descended  through  a  distance  vt/2, 
where  v  is  its  current  speed.  Thus,  an  equation  involving 
the  work  available  for  giving  rise  to  rotational  kinetic 


energy  of  the  plate  triangles  is. 


O  Vt/ .  i 

(p  -  p  )  4a^  ‘  - -  =  4  •  t  Iw2 

S  3  * 

Hence  substituting  for  p  from  (52) ,  noting  that  I 

S 


ma*/6 


and  simplifying, 
6M 

(p  -  -~§->  t  = 


/2  or  v 


(P-PJ  t 


or  directly  by  momentum  consideration,  thus, 


(  <p-ps  ) 


t 


a*  2a.  a  _  ma  v 
“1  ’  T  T  *  a 


(55) 


The  deflection  w^  at  the  end  of  time  t  is, 
wx  =  vt/2  =  (p  -  ps)T2/m 


(56) 


(2)  Medium  Load:  x  >  t 

After  the  removal  of  pressure  p  at  t  =  x,  further  deflec¬ 
tion,  w 2,  is  added  to  w^,  because  the  rotational  kinetic 
energy  of  the  four  plate  segments  must  be  dissipated  in  doing 
plastic  work  in  the  hinges.  Let  the  segments  rotate  through 
further  angle  <j>  before  coming  to  rest,  then, 

1  v2  w  / 

4  •  (=■  I  * — =>)  =  4  a/2  •  M  (<J>/2)  =  8aM  •  ^'a 

2  a2  P  P 


Hence 


W—  » 


maV 
2  4  M. 


™  2 
mv 

*57 


<P-PS>  2t2 
ra  Pe 


(57) 


using  Eqn  (55) .  Thus  the  total  central  deflection  is. 


w  -  wL  +  w2  =  p  (p  -  ps)x2/mps  (58) 

If  the  total  response  time  of  the  plate  is  xt,  note  that 
xtPg  *  px.  In  the  case  of  thin  plates  clamped  at  their  outer 
periphery  it  is  clear  that  after  a  significant  amount  of 
deflection  there  must  be  stretching.  An  analysis  of  the 
dishing  of  peripherally  clamped  metal  plates  when  given  a 
transverse  initial  speed  of  a  few  hundred  feet  per  second 
has  been  given  by  Hudson  [Ref.  3].  Descriptively,  Hudson's 


mu.  mw.uMWA’WW  '>.J .  y  rrr  r  -.  iv.- 


model  is  easily  understood  by  reference  to  Fig.  3-2.  Central 
portions  of  the  plate  maintain  their  given  transverse  speed 
until  overtaken  by  a  bending  wave  propagated  radially  inwards 
from  the  clamping  ring.  Fig.  3-2 (a)  shows  the  imagined  plate 
configuration  at  successive  intervals  during  the  deformation 
process;  when  the  plate  deflection  is  significant,  the  process 
is  predominantly  one  of  stretching  and  for  many  practical 
purposes  the  technological  approach  is  adequate. 

Figure  3-2 (a)  shows  the  terminal  state  taken  up 
by  a  thin  rectangular  steel  plate  clamped  around  its  periphery 
after  being  subjected  to  the  uniform  but  intense  impulsive 
loading,  the  maximum  speed  of  the  plate  is  a  few  hundred 
feet  per  second,  but  contracting  in  area  with  time;  the 
dynamic  situation  during  deformation,  with  plastic  hinges 
moving  into  the  plate  from  the  clamped  edges  and  intersecting 
along  lines  equally  inclined  to  the  sides  from  the  corners, 
is  indicated  diagrammatically  in  Fig.  3-2 (b) . 

Now  the  yield  criteria  of  metals  should  be  con¬ 
sidered.  Two  approaches  to  the  interpretation  of  a  yield 
criterion  are  possible;  one  is  Treses' s  criterion,  the  other 
is  Von  Mises  criterion.  But  for  most  metals  it  is  generally 
found  that  while  the  experimental  points  fall  between  the  two 
ellipses,  they  incline  towards  Von  Mises'  ellipse.  Experiments 
were  earlier  performed  by  Lode  (1925)  on  tubes  subject  to 
tension  and  internal  pressure.  Bending  and  torsion  was  used 
by  Siebel  in  1953.  All  results  confirm  the  general  opinion 
that  the  best  simple  yield  criterion  for  metals  is  that  of 
Von  Mises  [Ref.  5].  . 


a  ^  -V, 


B.  PROBLEM  DESCRIPTION  AND  NASTRAN  DECK  SETUP 
X.  Problem  Description 

The  model  is  shown  in  Fig.  3.3.  It  is  a  20  in.  by 
20  in.  steel  square  plate  model  which  is  divided  to  25  QUAD4 
(NASTRAN  card)  elements  and  its  boundary  is  clamped.  Its 
thickness  is  0.8  inches  and  the  engineering  data  as  follows. 

E  =  3.0E+7 
v  =  0.3 
L  =  20.0  in. 
p  =  7.33E-4  lb/in3 
t  =  0.8  in. 

where  E  is  the  elastic  modulus/  v  is  the  Poisson's  ratio, 

L  is  the  length  of  the  plate,  p  is  the  mass  density,  and 
t  is  a  thickness  of  the  plate. 

Next  the  shock  wave  applied  impinges  on  the  steel 
plate  in  the  -Z  direction  at  time  0.0  seconds  (Fig.  3.4)  and 
its  peak  overpressure  is  3.0E+4  psi  and  the  shock  wave  dura¬ 
tion  time  (t^)  is  5.0E-3  secs.  To  find  the  pressure,  the 
following  formula  is  used  [Ref.  6] 
p  =  Pm^^V  exP(“at/td) 

where  pm  is  peak  reflected  overpressure  at  time  0.0  secs 
and  a  is  wave  form  parameter.  But  the  shock  wave  pressure 
cannot  be  applied  as  a  load  to  the  plate  since  the  front  face 
shock  pressure  of  the  plate  is  relieved  by  the  rear  face  shock 
wave  which  is  positive  Z  direction.  The  following  formula  is 
used  [Ref .  4  ] . 


V  w>  .^v-  T 'X  '..*- 1  A',A’*‘V  ".*-’.,  ".<  ',.*^r.' 


2.  NASTRAN  Input  Data  Deck 

Reference  is  made  to  the  echo  of  the  input  data  deck 

which  appears  in  the  output. 

a.  Sol  66  [Ref.  6] 

1  NASTRAN  PREF0PT=2 

2  ID  NONLINEAR,  STATIC 

3  SOL  66 

4  TIME  30 

5  CEND 

6  TITLE-NONLINEAR  DEFORMATION  WITH  STATIC  LOADS 

7  SUBTITLE=TRANSIENT  RESPONSE  BY  SHOCK  WAVE  LOAD. 

8  ECHO=UNSORT 

10  DISP=ALL 

11  OLOAD=ALL 

12  SPCF=ALL 

13  GPFOR=ALL 

14  ELS T= ALL 

15  SEALL=ALL 

16  SPC=200 

17  SUBCASE  1 

18  LOAD=10 

19  NLPARM=10 

20  BEGIN  BULK 

21  $ 

22  $  LATERALLY  LOADED  FLAT  PLATE 

23  $ 

51 
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A-V 


24 

MESHOPT , , , , , YES 

25 

$ 

26 

$  USE  MSGMESH  TO  GENERATE  5*5  ARRAY  OF  QUAD4 

ELEMENTS 

27 

$ 

28 

EGRID, 1, ,0 . ,0  . 

29 

EGRID, 2 , ,20 . , 0 . 

30 

EGRID, 3, ,20. ,20. 

31 

EGRID, 4, ,0. ,20. 

32 

GRIDG, 2 , , , 5 , -1 , -2 , -3 , , +FLD2 

33 

+FLD2 ,5,-4 

34 

CGEN ,QUAD4 , 201,200 , 2 

35 

$  REFER  TO  PAGE  2 . 18-7 .APPLICATION  MANUAL. 

36 

PSHELL, 200 ,65,0.8,  65,, 65 

37 

MATl ,65,3. OE+7 ,,.3,7. 33E-4 

38 

MATS1, 65, 40, PLASTIC, ,1,1,2. 6E+4 

39 

TABLES  1,40,0, ,,,,, ,+TABl 

40 

+TAB1, 0 . 0 , 0 . 0 , 0 . 001 , 3 . OE+4 , 0 . 025 , 4 . -E+4 ,0.05 

,  5 . 8E+4 , +TAB2 

41 

+TAB2 ,0.10,6. 7E+4 ,0.125,7. OE+4 ,0.2,7. 5E+4 , FNDT 

42 

PARAM, COUP MASS , 1 

43 

PARAM , K6  ROT ,1.0 

44 

PARAM, LGDISP.l 

45 

PLOAD2 , 10 , -14 . 7 , 208 , 212 , 213 , 214 , 218 

46 

NLP ARM, 10 , 1, , AUTO , , 20 ,W, YES 

47 

$NLPARM, 20 , 3 , , AUTO, , 20 , W, YES 

48 

SPCG, 20 0,2, 123456, A, B 

49 

SPCG, 200, 2, 123456, B,C 
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card  28-31  -  Defines  vertices  of  quadrilateral  field 
using  MSGMESH  EGRID  command, 
card  33  -  Defines  a  field  of  grid  points  and  generate 
grid  cards. 

*  The  grid  number  generated  by  MSGMESH  starts 
from  the  FID  of  GRIDG  (if  FID  is  2,  the  first 
grid  point  ID  is  20000) 

*  5  means  5  elements  in  each  direction, 
card  34  -  Generate  element  connection. 

*  QCJAD4  is  the  element  type  which  should  be 
generated. 

*  Element  ID  number  starts  from  201. 

*  Fifth  field  200  is  the  PID  of  pshell  card, 
card  38  -  Specifies  table  references  and  material 

properties  which  are  stress-dependent  for  use 
in  material  nonlinearity  applications 

*  Taken  Von  Mises  yield  function,  isotropic 
hardening  rule. 

card  39  -  Represent  stress-strain  data. 

card  43  -  The  value  1.0  suppresses  singularities 

*  It  is  intended  primarily  for  geometric 
nonlinear  analysis. 

card  44  -  The  value  1  allows  the  nonlinear  calculations, 
card  45  -  Request  the  initial  condition  on  the  elements 
in  the  negative  Z  direction, 
card  46  -  Controls  the  nonlinear  analysis  iteration. 
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*  To  reduce  the  output,  use  sixth  field. 

*  For  rapid  convergence,  increase  the  Error 
tolerance  and  maximum  iteration  and  the  limit 
on  diverging  iterations. 

card  48-51  -  Generate  spcl  cards. 

This  card  is  only  used  to  MSGMESH  generation, 
b.  Sol  99 
NASTRAN  PREF0PT=2 
ID  NONLINEAR, TRANSIENT 
SOL  99 
TIME  59 
CEND 

TITLE=NONLINEAR  DEFORMATION  WITH  DYNAMIC  LOADS 
SUBTITLE=TRANSIENT  RESPONSE  BY  SHOCK  WAVE  LOAD. 
ECHO=SORT 

SET  1=20202,20203,20204,20304,20404 

SET  2=201,208,212,213,214,219,225 

SET  3=20000,20003,20005,20200,20400 

SDISP=1 

OLOAD=l 

SPCF=3 

ELST=2 

SEALL=ALL 

SPC=200 

LOADSET=77 

SUBCASE  1 
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X  . 


V 


20  DLOAD=  3  7 

21  TSTEPNL=45 

22  SUBCASE  2 

23  DLOAD=37 

24  TSTEPNL=46 

25  SUBCASE  3 

26  DLOAD=37 

27  TSTEPNL=47 

28  OUTPUT (XYPLOT) 

29  PLOTTER  NAST 

30  CSCALE  1 . 3 

31  XTITLE=TIME  IN  SECS 

32  XGRID  LINES=YES 

33  YGRID  LINES-YES 

34  YTI TLE= D I S PLACEMENT  GRID  20202 

35  XYPLOT  SDISP  RESP/20202 (T3) 

36  YTITLE** VELOCITY  GRID  20202 

37  XYPLOT  SVELO  RESP/20202 (T3) 

38  YTITLE-DISPLACEMENT  GRID  20101  20102  20202 

39  XYPLOT  SDISP  RESP/20101 (T3) , 20102 (T3) , 20202 (T3) 

40  BEGIN  BULK 

41  $ 

42  $  LATERALLY  LOADED  FLAT  PLATE 

43  $ 

44  MESHOPT, , , , ,YES 

45  $ 
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46  $  USE  MSGMESH  TO  GENERATE  5*5  ARRAY  OF  QUAD4  ELEMENTS 

47  $ 

48  EGRID,1,  ,0.  ,  0. 

49  EGRID,2, ,20. ,Q. 

50  EGRID,3, ,20. ,20. 

51  EGRID,4, ,0. ,20. 

52  GRIDG,2, , ,5, -1,-2, -3, ,+FLD2 

53  +FLD2 ,5,-4 

54  CGEN,QUAD4 ,201,200,2 

55  $  REFER  TO  PAGE  2 . 18-7 , APPLICATION  MANUAL 

56  PSHELL, 200, 65, 0.8,65, ,65 

57  MATl, 65 , 3 . OE+7 , , . 3, 7 . 33E-4 , , ,0 . 01 

58  MATS1, 65, 40, PLASTIC, ,1,1,2. 6E+4 

59  TABLESl, 40 , 0 ,,,,,, , +TAB1 

60  +TAB1 ,0.0,0.0,0.001,3. OE+4 ,0.025,4. OE+4 ,0.05,5. 8E+4 , +TAB2 

61  +TA32 ,0.10,6. 7E+4 ,0.125,7. OE+4 ,0.2,7. 5E+4 , ENDT 

62  PARAM , COUPMASS , 1 

63  $  PARAM , AUTOSPC , YES 

64  PARAM, W4,1.616E+3 

65  PARAM, K6 ROT, 1.0 

66  PARAM,  LGDISP , 1 

67  PARAM,  SLOOP I D, 1 

68  PARAM , LOOP ID , 1 

69  PARAM,  STIME, 0.000450 

70  PARAM ,  MAXLP ,  7 

71  DLOAD, 37, 1.0, 1.0, 7, 1.0, 6, 1.0, 5 
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TLO ADI, 7,61, ,0,22 
LSEQ, 77, 61,71 
PLOAD2, 71,3 .2E+4, 213 
TL0AD1,6,62, ,0,22 
LSEQ, 77, 62, 81 

PLOAD, 81, 2 . 5E+4 ,207,208,209,212,214 
TLO ADI, 5, 63, ,0,22 
LSEQ, 77, 63, 91 

PLOAD2, 91, 2 . 5E+4 ,217,218,219 
TABLED1 , 22 ,,,,,,, , +NEXT1 

+NEXT1,-1. ,0.0, 0.0, -0.0010, 0. 0001, -1.0, 0.0010, -0.560, +NEXT2 

+NEXT2, 0.004,-0. 001, ENDT 

TSTEPNL, 45 , 150 , 1. 0E-6 , 20 , AUTO, , 20 , W, +LL1 

+LL1,1. OE-7 , 1. OE-7 , 1. OE-7 ,20,20,10 

TSTEPNL, 46, 200, 1.0E-6, 20, AUTO, ,20, W,+LL2 

+LL2 ,1. OE-5 , 1. OE-5 , 1. OE-5, 20,20, 10 

TSTEPNL, 47, 120, 3. OE-6, 20 ,AUTO, ,20,W,+LL2 

+LL2 , 1. OE-4 ,1. OE-4 ,  1. 0E-4,20,20,10 

SPCG, 200 , 2 , 123456 , A,B 

SPCG, 200, 2, 123456, B,C 

SPCG, 200, 2, 123456, C,D 

SPCG, 200, 2, 123456, A, D 

ENDDATA 

1)  Executive  control  cards 

card  3  :  Requests  nonlinear  transient  analysis.  Rigid 


format  99 


2)  Case  control  cards 

card  18:  Selects  a  static  load  set  for  use  in  dynamics. 

Thus  this  may  be  referenced  by  dynamic  loads, 
card  20 :  This  card  selects  a  dynamic  load  for  nonlinear 
transient  problem. 

card  22-24:  Requests  restart  at  time  0.00015  secs, 
card  25-27:  Requests  restart  at  time  0.00035  secs, 
card  28-39:  Requests  the  xy  plotting  which  is 
corresponding  to  the  time  step. 

3)  Bulk  data  cards 

card  64:  This  gives  damping  effects  for  transient 
analysis . 

*  Sometimes,  more  rapid  convergence  effects 
can  be  obtained  by  this  damping  factor  with 
damping  coefficient  in  MATl  card  (GE) . 

card  67:  This  card  identifies  the  initial  conditions 
from  a  previous  SOL  66  nonlinear  static 
solution . 

*  Refer  to  the  par am  card,  SLOOPID. 

card  68,69:  Request  these  two  cards  when  we  taken  more 
than  one  subcase. 

card  70:  Requests  more  internal  calculation. 

*  We  can  reduce  iteration  steps  and  prevents 
the  overprinting  due  to  iteration  steps. 


card  71:  Defines  a  dynamic  loading  condition  for 

transient  response  problems  as  a  linear 

combination  of  load  sets 

*  See  the  remarks  5,  7,  9  of  DLOAD  card  in 
bulk  data. 

card  72,  75,  78:  Defines  time  dependent  dynamic  load. 

*  Requests  a  static  load  card  set  (LESQ) . 

card  73,  76,  79:  Defines  a  sequence  of  static  load  sets. 

This  load  sets  may  be  referenced  by  dynamic 
load  card  (TL0AD1) 

*  This  card  should  be  used  when  we  want  to 
apply  area  pressure  on  elements  than  a  point 
pressure  on  grid  point. 

card  74,  77,  80:  Defines  a  uniform  static  pressure  load 
on  the  plate  elements. 

card  81-83:  Defines  a  tabular  function  for  use  in 

generating  a  time-dependent  dynamic  loads. 

*  In  this  card  the  peak  pressure  (3.0E+4  psi) 
occurred  at  0.0001  sec. 

card  84,  86,  88:  This  card  control  the  nonlinear 
transient  analysis. 

*  This  card  is  expalined  in  the  sensitivity 
analysis . 

C.  SENSITIVITY  ANALYSIS 
1.  General 

In  transient  analysis,  two  types  of  instability  could 

occur.  The  first  is  the  familiar  nonlinear  load  iteration 

divergence  which  also  occurs  in  static  analysis.  The  second 
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is  the  divergence  which  grows  with  time  in  the  transient 
integration.  Both  instabilities  are  caused  by  uncorrected 
nonlinear  equilibrium  errors . 

A  convenient  method,  from  Von  Neumann,  for  analyzing 
the  stability  limits  is  to  assume  that  the  nonlinear  forces 
are  nearly  linear  and  the  error  vector  has  a  constant  conver¬ 
gence.  It  is  assumed  that  the  nonlinearity  has  a  first  order 
approximation:  [Ref.  7] 

(P(U)}  ~  [Kr  +  AKNL]  (u>  (59) 

where  AKNL  is  the  difference  between  the  tangent  stiffness 
matrix  and  its  approximation  Kr.  The  error  vectors,  {5}, 
are  assumed  to  grow  at  the  rate  X,  defined  as: 

{51+1}=  X{5i>  (60) 

Note  that  if  j  X  j  >  1  the  system  will  be  defined  as  unstable. 

After  lengthy  calculations  we  may  summarize  the  various 
criteria,  assuming  that  the  matrices  are  reduced  to  equivalent 
scalar  modal  quantities,  the  criterian  for  stable  solutions 
are: 

1)  For  time  steps  with  converged  "static"  iterations: 

1/2  >  0  >  1/4 

This  is  the  same  criteria  as  linear  analysis. 

2)  For  "static  load"  iterations: 

6  A  KNL  <  6  Kr  +  B  +  -iy  M 

This  states  that  the  mass  and  damping  add  to  the  effective 
linear  stiffness  [Kr]  and  improve  the  stability. 
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3)  For  time  step  integration  with  no  intermediate  static 
iterations,  as  in  the  standard  Newmark  Beta  method: 

A  KNL  <  (46  -1)  Kr  +  M 

At 

This  restriction  is  more  severe  than  criteria  b  above, 
proving  that  the  internal  iterations  are  more  stable 
than  the  Newmark  integration. 

In  summary,  the  Von  Neumann  method  will  have  fewer 
divergence  problems  than  either  the  static  nonlinear  solu¬ 
tion  or  the  single  step  transient  nonlinear  methods.  The 
better  stability  of  the  method  and  the  capability  to  use 
larger  time  steps  outweighs  the  cost  of  a  few  internal 
iterations  on  the  static  element  forces. 

2.  Nonlinear  Transient  Controls 

The  input  data  required  for  SOL  99  is  a  combination 
of  transient  control  data,  similar  to  SOL  69  (Direct  Transient, 
superelements) ,  and  nonlinear  modeling  data  similar  to  SOL  66 
(Nonlinear  Statics) .  The  nonlinear  properties  are  defined  by 
nonlinear  material  data  ( MATS 1 ) ,  gap  elements  (GAP) ,  and 
large  displacement  effects  (parameter  LGDISP) .  The  transient 
effects  are  produced  by  loading  functions  (TLOADi,  DAREA,  etc.) 
damping  (parameters,  elements,  and  material  data) ,  and  element 
mass  properties. 

The  only  unique  data  required  for  SOL  99  is  supplied 
on  the  TSTEPNL  data  cards.  (Several  optional  parameters  are 
provided  for  tuning  the  method  and  for  use  in  restarts.)  The 
TSTEPNL  card  in  itself  is  an  amalgam  of  the  normal  transient 
TSTEP  card  and  the  nonlinear  control  card,  NLP ARM. 


a.  Case  Control 

The  significant  rules  for  nonlinear  transient  are: 

1)  Each  subcase  defines  a  time  interval  starting  from  the 
last  time  step  of  the  previous  subcase  which  is  subdivided 
into  small  time  steps.  The  output  TIME  printout  will  be 
labeled  by  the  net  time,  including  all  previous  subcases. 

2)  The  input  loading  functions  may  be  changed  for  each  sub¬ 
case  or  continued  by  repeating  the  same  DLOAD  request. 

The  bulk  functions  are  evaluated  at  the  net  model  time, 
which  starts  at  the  last  time  step  of  the  previous  subcase. 

3)  Each  subcase  may  have  a  different  time  step  size,  time 
interval,  and  iteration  control  selected  by  the  TSTEPNL 
request.  The  case  control  requests  which  may  not  be 
changed  after  the  first  subcase  are:  SPC,  MPC ,  DMIG,  and 
TF. 

4)  Output  requests  for  each  subcase  are  processed  indepen¬ 
dently  and  the  output  for  each  case  are  generated  before 
the  next  subcase  is  processed.  Separate  XYPLOT  ouputs  will 
occur  for  each  subcase.  Available  outputs  are  DISPLACEMENT, 
VELOCITY,  ACCELERATION,  OLOAD,  STRESS,  FORCE,  SDISPLACEMENT, 
SVELOCITY,  and,  SACCELE RATION .  Note  that  NONLINEAR  and 
SPCFORCE  outputs  are  not  provided. 

5)  Initial  conditions  are  not  explicitly  provided  for  user 
input  in  SOL  99.  The  reason  is  that  the  nonlinear  element 
configuration  and  equilibrium  state  must  be  known  at  every 
time  step.  If  initial  conditions  were  given,  all  of  the 
nonlinear  element  forces  and  plastic  stresses,  as  well  as 
the  past  history,  must  be  given.  Therefore,  the  initial 
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conditions  become  part  of  the  analysis  and  are  generated 
by  initial  transient  subcase  or  by  restarts  from  a  previous 
static  solution. 

6}  All  normal  superelement  model  generation  options  and 

matrix  reduction  options  are  allowed  for  the  linear  portion 
of  the  structure.  General  Dynamic  Reduction,  Component 
Mode  Synthesis,  and  Guyan  Reduction  may  be  performed  for 
nonresidual  super elements.  Therefore,  the  residual  super¬ 
element  may  contain  scalar  degrees  of  freedom  representing 
linear  model  formulations . 
b.  TSTEPNL  Data 

The  TSTEPNL  data  card  description  is  provided  in 
Section  2  of  the  MS C/NAS TRAN  User's  Manual.  The  input  fields 
specify  the  time  step  size,  the  number  of  steps,  and  the  out¬ 
put  interval  as  well  as  the  nonlinear  iteration  options. 

The  nonlinear  iteration  data  are  provided  for 
controlling  the  accuracy  and  stability  of  the  nonlinear 
calculations  and  they  will  not  effect  linear  solutions.  The 
controls  are  similar  to  those  on  the  NLP ARM  bulk  data  used  in 
static  nonlinear  analysis. 

As  discussed  above  two  types  of  errors  or  insta¬ 
bility  are  encountered  in  the  transient  solution.  The  load 
iterations  are  similar  to  static  load  analysis  and  may  diverge 
at  one  or  more  time  steps  when  large  nonlinear  stiffness 
changes  occur.  Load  iteration  divergence  may  be  cured  by 
smaller  time  steps  or  by  stiffness  matrix  updates.  The 


50)  indicates  the  rate  of  growth  of  this  error. 


The  size  of  the  time  step  is  of  primary  importance 


PARAM,SLOPPID,N 


where  N  was  the  printed  value  of  LOOP  for  the  desired 
static  solution.  Constraint  sets,  direct  input  matrices, 
mass,  and  damping  may  be  changed. 

Restarts  from  a  previous  transient  execution  are 
available  for  a  number  of  cases.  If  the  same  model  is  to  be 
reexecuted,  only  the  residual  superelement  needs  to  be  re¬ 
assembled  (SEMA,  SELA=0 ) .  If  the  end  results  for  the  previous 
transient  run  are  to  be  used  for  the  initial  conditions  at 
t=o,  add  the  parameter  L00PID=N  where  N  corresponds  to  the 
Nth  subcase  of  the  previous  run  and  N  dummy  SUBCASE  commands 
to  start  the  residual  case  control  execution. 

The  normal  restart  for  a  transient  run  is  to 
continue  from  the  last  step  of  a  previous  subcase  with 
different  loads  and/or  TSTEPNL  data.  Provide  the  following 
parameters : 

LOOPID  =  N  -  Start  from  the  Nth  subcase 

Stime  =  t  -  Start  from  time  t. 

Note  that  constraint  sets  should  not  be  changed 
to  avoid  incompatible  matrix  sizes. 

D.  NUMERICAL  RESULTS 

The  progressive  deformations  of  the  plate  are  plotted 
along  the  x-axis  at  five  different  time  locations  as  shown 
in  Fig.  3.6(b).  At  the  early  time  (1.5E-4  secs),  the  response 
of  the  plate  is  in  the  elastic  region.  During  the  deformation, 
plastic  hinge  lines  are  formed  along  the  clamped  edges  and  the 
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diagonals  of  the  plate  as  indicated  in  Fig.  3.6(a).  This 
response  phenomenon  agrees  well  with  Hudson's  studies 
[Ref.  4].  As  described  in  Section  A-4  of  Chapter  3,  when 
the  plate  deflection  is  significant,  the  plate  experiences 
in-plane  stretching. 

Time  history  displacement  responses  of  Grid  points 
20101,  20201,  and  20202  are  plotted  in  Fig.  3.7.  As  shown 
in  the  figure,  for  the  short  duration  impulse,  the  maximum 
displacement  occurs  at  the  center  of  the  plate  after  termi¬ 
nation  of  shock  loading.  The  rate  of  increase  in  displacement 
is  the  maximum  at  the  termination  of  shock  loading  as  shown 
in  Fig.  3.8.  During  the  application  of  the  shock  loading, 
the  acceleration  level  is  high,  implying  the  high  inertia 
force  in  the  system.  At  the  termination  of  the  inertia 
force  in  the  system,  the  acceleration  stays  in  the  constant 
level  until  the  plate  displacement  reaches  the  maximum. 

Time  h iso try  Von-Mises  stress  of  QUAD4  element  213  (plate 
center)  is  plotted  in  Fig.  3.7.  The  stress  variation  is 
highly  non-linear  due  to  the  material  non-linearity. 


Figure  3.6(b).  Progressive  Deformation  of  the  Plate 

Along  Section  A-A. 


Figure  3.7.  Time  History  Displacement  Responses  of  Grid  Points 
20101,  20201,  20202. 


Figure  3.8.  Time  History  Velocity  Responses  of  Grid  Point  20202 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  state-of-the-art  literature  review  revealed  that  the 
skirt  plate  is  very  effective  to  absorb  initial  explosion 
energy  before  the  shock  loading  reaches  the  main  armor.  The 
material  of  the  armor  must  possess  the  high  hardness  to  prevent 
the  severe  damage  and  the  high  ductility  to  minimize  the 
spalling.  The  composite  materials  may  serve  more  effectively 
them  the  conventional  material  to  minimize  the  damages  against 
the  powerful  modern  projectiles. 

The  non-linear  finite  element  analysis  of  the  plate  sub¬ 
jected  to  the  air  shock  loading  verifies  the  plastic  hinge 
line  formation  and  in-plane  stretching  inducing  snap-through 
collapse.  For  the  short  duration  shock  loading,  the  maximum 
deformation  occurs  after  the  termination  of  loading,  and  the 
velocity  reaches  the  maximum  at  the  termination  of  shock. 

The  Von-Mises  stress  is  highly  non-linear  due  to  the  material 
non-linearity.  For  the  non-linear  analysis  using  MSC/NASTRAN , 
TSTEPNL  bulk  data  card  is  found  to  be  the  most  important  card 
which  controls  the  accuracy  and  the  stability  of  non-linear 
calculation. 

To  insure  the  convergence  of  solutions,  convergence  limit 
in  the  TSTEPNL  card  in  the  bulk  data  deck  should  be  considered. 
If  a  program  was  run  with  big  convergence  limit,  the  solution 
will  be  gotten  rapidly  with  small  iteration  steps  but  it  is 


not  accurate  and  large  time  steps  can't  be  taken.  On  the 
contrary,  if  the  program  was  run  with  small  convergence 
limit,  the  solution  will  be  accurate  but  it  is  gotten  slowly 
with  large  iteration  steps.  Thus,  it  is  expensive  and  anyway 
large  time  steps  can't  be  taken  too.  Thus,  compromise  should 
be  taken  to  run  a  large  time  steps.  For  example,  if  one  non¬ 
linear  transient  problem  is  needed  total  700  time  steps  to  get 
the  solutions,  the  following  procedure  is  recommended.  In  the 
first  subcase,  the  program  is  iterated  with  smaller  conver¬ 
gence  limit,  10*E-7  to  get  more  accurate  solution  and  in  the 
second  subcase,  the  program  is  run  with  limit,  lO*E-5,  and  in 
the  third  subcase,  the  program  is  iterated  with  convergence 
limit,  10*E-4. 

Consequently,  the  total  700  time  step  will  be  calculated 
successively. 

The  sensitivity  studies  on  shock  period  and  amplitude  is 
suggested  to  understand  their  effects  on  plate  responses. 


APPENDIX  A 


THE  COMMAND  TO  OPERATE  DATA  BASE 

****THE  COMMAND  TO  ERASE  THE  DATA  BASE  NAMED  "MSS .S2770P .NODE" 

//ERASE  JOB  (2770,0252) ,NASTRAN,CLASS=A 
//♦MAIN  ORG=NPGVMl . 2 7 70P 
//  EXEC  PGM=IEFER14 
//SYSPRINT  DD  SYSOUT=A 

//DDl  DD  DSN=MSS.S2770. NODE, DISP=(OLD, DELETE) 

// 

****THE  COMMAND  TO  FLOT  FROM  THE  DATABASE  NAMED  "MSS .S2770 .PPLOT" 

//LEEPLCT  JOB  (2770 , 0252) ,NASTRAN,CLASS=C 
//♦MAIN  ORG=NPGVM1.2770P 
//♦FORMAT  PR,DDNAME=SYSVECTR, DEST=LOCAL 
//  EXEC  NASTPLOT,PLTDSN='MSS.S2770. PPLOT' 
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APPENDIX  B 


DATA  DECK  FOR  PROGRAM 


GENERAL  DESCRIPTION  Of  DATA  OEOC 

The  Input  deck  begins  with  the  required  resident  operating  system  control  cards.  The  type 
and  nutter  of  these  cards  will  vary  with  the  installation.  Instructions  for  the  preparation  of 
these  control  cards  are  given  in  Section  7.6  of  the  MSC/NASTRAN  Application  Manual. 

The  operating  system  control  cards  are  followed  by  the  MSC/NASTRAN  Oata  Deck,  which  consists 
of  the  following  three  sections: 

1.  Executive  Control  Oeck 

2.  Case  Control  Oeck 

3.  Bulk  Oata  Deck 

The  Executive  Control  Oeck  and  the  Case  Control  Oeck  both  have  free- field  formats.  Only 
coliaans  1  through  72  are  used  for  data.  Any  Information  in  colons  73  through  80  may  appear  in 
the  printed  echo,  but  the  data  will  not  be  used  by  the  program.  As  explained  In  Section  2.4.1, 
limited  use  is  made  of  the  data  in  coliaans  73  through  80  for  the  Bulk  Oata  Oeck.  If  the  last 
character  on  a  card  is  a  conea  (not  necessarily  in  column  72),  the  next  card  Is  a  continuation  of 
this  physical  card.  Any  matter  of  continuation  cards  may  be  specified,  and  together  they  form  a 
logical  card. 

The  NASTRAN  card  is  used  to  change  the  default  values  for  certain  operational  parameters, 
such  as  buffer  size  or  the  number  of  data  lines  printed  per  page.  More  than  one  NASTRAN  card  may 
be  present.  NASTRAN  cards  are  optional,  but.  tf  present,  they  mist  be  the  first  cards  of  the  data 
deck.  The  NASTRAN  card  is  a  free- field  card  (similar  to  cards  in  the  Executive  Control  Oeck). 
Its  format  is  as  follows: 


NASTRAN  keyword^  ■  value,  keyword;  •  value,  ... 
An  alternate  format  is: 

NASTRAN  keyword}  •  value 
NASTRAN  keyword;  ■  value 
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Input  Data  Card  EGRIO  Grid  Point  Position 

Description:  Defines  the  position  of  a  geometric  grid  point  of  the  structural  model 


format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

3 

9 

10 

EGRIO 

ID 

a> 

XI 

X2 

X3 

X 

X 

X 

EGRIO 

2 

3 

1.0 

2.0  I 

3.0 

Field 

Contents 

10 

Grid  point  Identification  number  (Integer  >  0) 

CP 

Identification  number  of  coordinate  system  in  which  the 
Is  defined  (Integer  ^0  or  blank! 

position  of  the  grid  point 

X1.X2.X3 

Position  of  the  grid  point  In -coordinate  system  CP  (Real 

or  Blank) 

Remarks:  1. 

The  meaning  of  XI,  X2,  and  X3  depends  on  the  type  of 
follows:  (see  CBRD1  card  descriptions) 

coordinate  system,  CP,  as 

Type 

XI 

X2 

X3 

Rectangular 

X 

Y 

Z 

Cylindrical 

R 

9 (degrees) 

Z 

Spherical 

R 

9 (degrees) 

♦(degrees) 

2.  The  EGRIO  card  Is  used  In  MSGMESH  to  define  curved  edges  and  for  other  special 
purposes  (see  Part  I,  Section  3.0).  It  is  similar  to  a  GRID  card,  except  that  it 
will  not  be  processed  by  MSC/NASTRAN. 


3.  Entries  In  fields  7,  3  and  9  are  not  checked  by  NSGMESH. 


Input  Oats  Card  SHIPS  Linear  and  Higher  Order  Interpolation 


Description:  Defines  a  field  of  grid  points  and  generates  SKID  cards 
Format  and  Example: 


1  23456789  10 


GRIOG 

FIO 

CO/CI 

PS 

L 

GA 

GB  GC 

SEIO 

IQ . 

*0 

18 

-1 

-2 

-3 

*61 

M 

GO 

N 

GE 

GF 

GG 

GH 

Cl 

♦Gl 

9 

-4 

•H52 

El  (AS) 

E2(A8) 

a  (so 

£2(80 

El  ( C0 1 

E2(C0) 

El(OA) 

E2(DA) 

1 

♦G2 

6 

5 

8 

7 

El(EF) 

£2(EF) 

a(FG) 

E21FG) 

El(GH) 

E2(GH) 

El  ( HE) 

E2(  HE) 

El(AE) 

E1(8F) 

E2(  8F) 

El(CG) 

E2ICG) 

El(OH) 

E2(DH) 

Field 

FIO 

CO/CI 

PS 

L,M,M 

SEIO 

6A . SH 

Cl 

El(XT) 


Contents 

10  number  of  a  GRIDS  field  (Integer,  1  FIO  <_  99).  Must  be  unique  for  all  GRIDG 
cards.  When  entry  Is  negative.  Internal  EGRIO  cards  are  generated. 

Prescribes  displacement  and/or  Interpolation  coordinate  systems  (Integer  or 
blank).  See  Remarks  4  and  5. 

PS  entry  on  generated  GRID  cards  (any  of  the  digits  1-6  with  no  Imbedded  blanks) 
(Integer  ^0  or  blank) 

Grid  field  dimensions  when  an  entry,  L,  M,  and/or  N  Is  positive.  (See  Figure  3.) 
Absolute  value  is  a  LID  entry  on  a  LIST  Input  card,  when  an  entry  Is  negative. 

SEIO  entry  on  generated  GRID  cards  (Integer  >_0  or  blank) 

Grid  point  Identification  numbers,  GIDs,  of  grid  field  vertex  points  (Integer  A  0  or 
blank).  See  Remarks  3  and  6. 

CP  entry  on  generated  GRID  cards'  (Integer  ^  0  or  blank),  and  grid  field  inter¬ 
polation  coordinate  system  definition.  See  Remarks  4  and  5. 

10  number  of  edge  grid  point  (closer  to  X)  on  side  XT  (Integer  A  0  or  blank).  See 
Remarks  3  and  6. 


E2(XY)  ID  number  of  edge  grid  point  (closer  to  Y)  on  side  XY  (Integer  A  0  or  blank).  See 

Remarks  3  and  6. 


(Continued) 


GRIDS  (Cont.) 


Remarks:  1. 


Sec  Part  I,  Sections  2  and  3  for  discussions  on  the  GRIDS  Input  card. 


2.  Required  data  by  type  of  field 


Type 

Required  Data 

LIKE 

FIO.  L,  GA,  SB 

TRIA 

■  FID,  l.  GA,  G8,  GC 

QUAD 

FID,  L.  H,  GA.  GB,  GC.  GO 

HEX 

FID,  L,  M,  N,  GA  thru  GH 

Only  the  parent  and  first  continuation  card  are  applicable  for  linear 
Interpolation. 

3.  Vertex  and  edge  points  are  defined  on  (SR  10  or  EGRID  Input  cards,  or  on  GRID  cards 
generated  by  previously  processed  GRIDS  Input  cards.  Field  vertex  points  are 
automatically  equlvalenced  to  corresponding  entries  GA,  ....  GH,  when  these  entries 
are  positive.  GRID  cards  are  generated  only  for  active  field  grid  points.  See 
Part  I,  Section  10. 1. 

4.  CO /Cl  entry  may  be  used  to  prescribe  both  the  displacement  coordinate  system,  CD, 
and  the  Interpolation  coordinate  system.  Cl.  Both  CO  and  Cl  are  prescribed  when 
the  entry  Is  positively  signed.  Only  CD  Is  prescribed  when  the  entry  Is 
unsigned.  Only  Cl  Is  prescribed  when  the  entry  Is  negatively  signed.  The  Cl  entry 
on  the  first  continuation  card  may  be  used  when  both  CO  and  Cl  exist  and  are 
dl f ferent. 

5.  Cl  references  a  CIO  entry  on  a  coordinate  system  definition  Input  card.  When  Cl  Is 
entered  as  0,  Interpolation  is  In  basic  rectangular  coordinates.  When  Cl  Is  not 
entered.  Interpolation  Is  in  the  coordinate  system  that  positions  the  field  vertex 
and  edge  points.  In  this  case,  all  vertex  and  edge  points  mist  b*  defined  In  the 
same  coordinate  system  (CP  entries  on  GRID  and/or  EGRID  cards).  See  Part  I, 
Section  3.1  and  Table  2. 

6.  The  a  Interpolation  point  components.  In  spherical  coordinate  systems,  or 

the  8  Interpolation  point  components.  In  tyllndrlcal  coordinate  systems,  may  be 
offset  by  a  t360*  by  appending  a  tS  or  a  tC  to  grid  point  entries. 

7.  If  an  entry  for  GO  Is  present  In  a  TRIA  field,  grid  point  D  replaces  grid  point  C 

In  Interpolating  grid  point  positions  from  side  CA.  Grid  point  0  is  usually 
defined  by  an  EGRID  card.  N  must  be  blank.  See  Part  I,  Section  3.S. 

8.  If  both  El(XY)  and  E2 (XT )  are  entered,  the  positions  of  generated  grid  points  are 

determined  by  cubic  Interpolation  from  edge  XY.  If  E2(XY)  Is  blank,  positions  are 
determined  by  quadratic  Interpolation.  If  both  El (XY )  and  E2(XY)  are  blank, 

positions  are  determined  by  linear  Interpolation.  Either  none  or  one  edge  point, 

E(XY)  may  be  specified  for  each  edge  of  a  TRIA  field.  El(CD)  positions  a  grid 

point  on  edge  CA  of  TRIA  fields.  See  Part  I,  Section  3.7. 

9.  Intermediate  continuation  cards  cannot  be  omitted  even  If  all  entries  are  blank.  - 
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Input  Data  Card  CGEN  Generate  Element  Connection  Bulk  Data  Cards,  Form  1 


Description:  Generate  element  connection  cards  for  element  types  listed  under  Remark  1 
Format  and  Example: 


12  34S6789  10 


CGEN 

TYPE 

FEID 

PID 

FIO 

OIR 

TH 

. 

EIOL 

EIDH 

CGEN 

QUA04 

100 

10 

1 

♦<2 

.  IK' 

■TC-- 

T5 

X 

x 

xx 

+C2 

0.002 

0.002 

0.001 

0.001 

Field  Contents 

TYPE  Any  of  the  mnemonics  listed  under  Remark  1 


FEID 

PID 

FID 

DIR 

TH 


EIDL 


The  absolute  value  of  FEID  fs  the  element  Identification  number,  EID,  of  first 
element  In  element  set  (Integer  f  01 .  See  Remark  4. 

ID  of  an  element  property  card  or  PGEN  card  (Integer  >  0  or  blank).  Must  be  unique 
with  resoect  to  all  other  property  card  IDs  and  PGEN  IDs.  See  Remark  5. 

10  of  associated  grid  point  field  defined  on  a  GRIDG  or  GRIDU  card 

One  of  the  entries  l.  M,  or  H  (Default  entry  Is  L).  Direction  of  elements  In  the 
field.  See  Table  4  of  Part  I,  Section  7.1. 

Material  prooerty  orientation  specification  (Real  or  Blank;  or  Integer  >_  0  and  less 
than  1,000,000).  If  Real  or  Blank,  specifies  the  material  property  orientation 
angle  In  degrees.  If  integer,  the  orientation  of  the  material  x-axIs  Is  along  the 
projection  onto  the  plane  of  the  element  of  the  x-axIs  of  the  co-ordinate  system 
specified  by  the  Integer  value. 

ID  of  element,  below  which  no  element  connection  cards  are  generated  (Integer  >  0 
or  blank) 


EIDH 


ID  of  element,  above  which  no  element  connection  cards  are  generated  (Integer  >  0 
or  blank) 


TA.TB.TC.TO  Thickness  of  comer  points  of  qUAD  and  TRIA  fields  when  TYPE  Is  TRIA3,  TRIA6,  QUAD* 
or  QUA08  (Real  or  Blank).  See  Remark  6. 


Remarks:  1.  Element  types:  R0O,  P10TEL,  CFTUBE,  TUBE,  TRIA3,  TRAPRG,  TRIARG.  TRIAS,  TRIMS, 
TRIAX6,  QUA04,  SHEAR,  QUAD8,  HEXA1,  HEXA2,  HEXA8,  HEXA20  and  HEX20.  The  mnemonics 
HEXA8  and  HEXA20  refer  to  HDIA  elements  with  8  or  9  through  20  nodes,  respectively. 


(Continued) 
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,  O  -  '  •  *  .  • 


CGEN  (Cont.) 


2.  Elements  are  tabulated  by  class  In  Table  3,  of  Part  I,  Section  7.0.  Line  class 

elements  may  be  generated  In  all  field  types.  Tr1a3  and  Triad  class  elements  may 

be  generated  In  TRIA  and  QUAO  fields.  Quad4  and  Quad8  class  elements  may  be 
generated  In  TRIA,  QUAO  and  HOC  fields.  Only  QUA04  and  QUAD6  element  types  are 

valid  In  TRIA  fields.  TRIA3  and  TRIA6  elements  are  generated.  In  this  case,  along 

the  field  edge  correspond! ng  to  the  element  direction.  Hex8  and  Hex20  class 
elements  may  only  be  generated  In  HEX  fields. 

3.  MSGMESH  automatically  generates  PENTA  elements  In  place  of  HEXA  and  HEX20  elements, 
when  adjacent  element  edges  would  have  the  sane  grid  point  Identifiers.  TRIA6 
elements  replace  QUADS  elements  and  TRIA3  elements  replace  QUAD4  elements,  when 
adjacent  element  vertices  would  have  the  same  grid  point  Identifiers.  Otherwise, 
MSGMESH  suppresses  the  generation  of  an  element  which  would  reference  Identical 
grid  point  Identifiers.  Elements  are  generated  only  when  all  vertex  points 
exist.  Oeleted  midside  nodes  suppress  the  generation  fo  TRIM6  and  TRIAX6  elements. 

4.  EIDs  are  Incremented  when  the  FEID  entry  is  positive  and  decremented  when  the  FEID 
entry  Is  negative. 

5.  The  PID  entry  may  be  blank  for  those  elements  whose  element  connection  card  permit 
a  blank  property  10. 

6.  Element  thicknesses  may  be  generated  for  TRIA3,  TRIA6,  QUAD4  and  QUA08  elements 
located  In  TRIA  and  QUAO  fields.  Interpolation  from  the  thicknesses  at  the  field 
vertices  Is  based  on  field  topologies  with  equal  spaclngs  of  thicknesses  for  both 
GRIOG  and  GRIQ.U  fields.  Unequal  spaclngs  of  thicknesses  are  defined  for  GRIDG 
fields  specified  with  unequal  spaclngs  of  grid  points. 


Input  Oata  Card  MAT!  Material  Property  Definition,  Form  1 

Description:  Defines  the  material  properties  for  linear,  temperature-independent,  isotropic 

materials 


Format  and  Example: 
1  2 


MAT1 

MID 

MAH 

17 

ST 

*8C  20.+4 


Field  Contents 

MID  Material  Identification  number  (Integer  >  0) 

E  Young's  modulus  (Real  or  blank) 

G  Shear  modulus  (Real  or  blank) 

NU  Poisson's  ratio  (-1.0  <  Real  _<  0.5  or  blank) 

RH0  Mass  density  (Real) 

A  Thermal  expansion  coefficient  (Real) 

TREF  Thermal  expansion  reference  temperature  (Real) 

GE  Structural  element  damoing  coefficient  (Real) 

ST.SC.SS  Stress  limits  for  tension,  compression,  and  shear  (Real).  (Used  only  to  coroute 

margins  of  safety  In  certain  elements;  they  have  no  effect  on  the  computational 
procedures. ) 

MCSIO  Material  Coordinate  System  identification  number  (Integer  _>  0  or  blank) 


Remarks:  1.  The  material  Identification  number  must  be  unique  for  all  MAT1,  MATZ,  MAT3  and  MAr9 
^  cards. 

2.  MAT1  materials  may  be  made  temperature  dependent  by  use  of  the  MATTI  card. 

3.  The  mass  density,  RHO,  will  be  used  to  automatically  compute  mass  for  all 
structural  elements. 

4.  Weight  density  may  be  used  In  field  5  if  the  value  l/g  is  entered  on  the  param  card 
WTMASS,  where  g  Is  the  acceleration  of  gravity  (see  Section  3.1.5). 

5.  MCSIO  must  be  nonzero  If  the  CURV  module  Is  used  to  calculate  stresses  or  strains 
at  grid  points  on  plate  and  shell  elements  only. 

6.  To  obtain  the  damping  coefficient,  GE,  nultiply  the  critical  damping  ratio  C/CQ.  by 

2.0. 

(Continued) 
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MATl  (Cont.) 


7.  Either  E  or  G  must  be  specified  (l.e.,  nonolank). 

3.  If  any  one  of  E,  G,  or  NU  is  blank,  it  <*111  be  computed  to  satisfy  the  identity 

E  ■  2(1*NU)G;  otherwise,  values  supplied  by  the  user  will  be  used.  this 
calculation  Is  only  made  for  Initial  values  of  E,  G,  and  NU. 


9.  If  E  and  NU  or  G  and  NU  are  both  blank,  they  will  both  be  given  the  value  0.0. 

10.  Implausible  data  on  one  or  more  NATl  cards  will  result  In  a  warning  message. 
Implausible  data  is  defined  as  any  of  E  <  0.0,  G  <  0.0,  0.5  <  NU  <  0.0,  or 

|1  -  £ i I  >  0-01  (except  for  cases  covered  by  Remark  9). 


Input  Oata  Card  MATS1  Material  Stress  Dependence 


Description:  Specifies  table  references  and  material  properties  which  are  stress-dependent  for 
use  in  material  nor  linearity  applications.  This  card  will  be  activated  If  MATl  with  the  same  MID 
Is  being  used  In  the  material  nonlinear  solution  sequence  (66). 


Format  and  Example; 


12  3  4  5  6  7  a  9  10 


MATS1 

MID 

Rl 

TYPE 

3 

YF 

HR 

LIMITl 

x 

MATSI 

17 

28 

PLASTIC 

2 

3 

2. +4 

♦ABC 

LIMITZ 

♦ABC 

Field  Contents 

MID  Identification  number  of  a  MAT1  card  (Integer  >  0) 

R1  Identification  number  of  a  TA8LES1  card  ( Integer  ^  0) 

TYPE  Type  of  material  nonlinearity  (BCD)  NLELAST  (Nonlinear  elastic)  or  PLASTIC 

(Plasticity  theory) 

9  Wort  hardening  slope  (slope  of  stress  vs.  plastic  strain)  In  units  of  stress 

(Real).  For  elastlc-perfectty  plastic  cases,  B  »  0.0.  For  more  than  a  single 
slope  In  the  plastic  range,  the  stress-strain  data  must  be  supplied  In  a  TABLESI 
card  referenced  In  field  3.  See  Remarks  I  and  2. 

Yield  function  (Integer);  one  of  the  following  values:  (Default  *  I) 

1  (von  Miles) 

2  (Tresca) 

3  (Mohr-Coulomb) 

4  (Orucker-Prager) 

Hardening  Rule  (Integer);  one  of  the  following  values:  (Default  •  I) 

1  (Isotropic) 

2  (Kinematic) 

3  (Combined  Isotropic  and  kinematic  hardening) 


Parameter  used  In  the  yield  function  specification  (Real).  See  Remark  4. 

Parameter  used  In  the  yield  function  specification  for  yield  functions  3  and  4 
(Real).  See  Remark  4. 


MATS1  (Cont.) 


Remarks:  1.  rf  type  «  NLELAST ,  the  stress-strain  data  given  In  the  TA8LES1  card  will  be  use< 
determine  the  stress  for  a  given  value  of  strain.  The  values  8,  YF,  HR,  LIM 
and  IIMIT2  will  not  be  used  In  this  case. 

2.  If  type  *  PLASTIC,  either  the  table  Identification,  Rl,  or  the  work  hardei 
slope,  8,  may  be  given  but  not  both.  If  the  table  10  Is  omitted,  the  > 
hardening  slope,  8,  specified  In  Meld  S  Is  defined  as 


3 


or 


where  E  Is  the  elastic  modulus  and  f*  Is  the  slope  of  the  uniaxial  stress-str 
curve  In  the  plastic  region. 

3.  If  Rl  Is  given,  TA8LES1  entries  of  stress-strain  data  must  conform  to  the  follow 
rules: 


a.  If  type  ■  NLELAST,  the  curve  must  pass  through  origin  (x.-0.,  y,»0.,  for  s 
value  of  1).  This  option  Is  valid  only  for  line  elements. 

b.  If  type  ■  PLASTIC,  only  data  for  the  first  quadrant  need  be  supplied.  Th¬ 
ere  two  options  available  to  specify  the  data.  In  case  I  (defined  as  F0RI 
In  the  TABLES1  card)  the  first  point  Is  at  the  origin,  the  second  point  Is 
the  yield  stress  given  In  field  3  of  the  MATS1  card.  The  slope  of  the  1 
joining  the  origin  to  the  yield  stress  must  be  equal  to  the  value  of  E  g1« 
on  the  MAT1  Bulk  Data  card.  In  case  2  (defined  as  FORM- l  In  the  TABLI 
card)  the  slopes  of  the  curve  In  the  plastic  range  are  supplied.  See  t 
TABLES1  card  description  for  additional  details. 

4.  LIMIT!  and  LIMIT2  are  parameters  used  in  the  yield  function  definition  as  follows 


Yield  Function 

LIMIT! 

LIMIT2 

von  Mlses  (1) 

Yield  stress  In  tension,  <» 

Not  used 

or  Tresca  (2) 

y 

Mohr-Coulomb  (3) 

Cohesion, c  (In  stress  units) 

Angle  of  Internal 

or  Orucfcer-Prager  (4) 

friction  a  (In  degree: 

S.  See  Section  2.14.1.2  In  the  MSC/NASTRAN  Application  Manual  for  a  detail 
description  of  yield  function  and  hardening  rules. 
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Input  Data  Card  TABLES1  Material  Property  Table,  Form  I 

Oescrl ptlon;  Oeflnes  a  tabular  function  for  stress-dependent  material  properties  sucn  as  the 
stress-strain  curve  and  creep  parameters 


Field  Contents 

10  Table  Identification  number  (Integer  >  0). 

I  Not  used 

x^  ,y^  Tabular  entries  (Real) 


Remarks :  1.  The  x^  must  be  in  either  ascending  or  descending  order  but  not  botn. 

2.  Jumps  between  two  points  (x^  *  are  allowed,  but  not  at  the  end  points. 

3.  At  least  two  entries  must  be  present. 


4.  Any  x-y  entry  may  be  Ignored  by  placing  the  SCO  string  SKIP  In  either  of  the  two 

fields  used  for  that  entry. 

5.  The  end  of  the  table  Is  Indicated  by  tne  existence  of  the  BCD  string  ENOT  in  either 

of  the  two  fields  following  the  last  entry.  An  error  is  detected  if  any 

continuation  cards  follow  the  card  containing  the  end-of-table  flag  ENOT. 

6.  TA8LES1  Is  used  to  input  a  curve  in  the  form  of 


wnere  X  is  input  to  the  table  and  Y  is  returned.  The  table  look-up  yx(*).  x  ■  X, 
Is  performed  using  linear  Interpolation  within  the  table  and  linear  extrapolation 
outside  the  table  using  the  last  two  end  points  at  the  appropriate  table  end.  At 
4uap  points  the  average  yy(x)  Is  used. 


36 


Input  Oata  Card  PlflA02  Pressure  Load  on  a  Two-0 f mens Iona I  Structural  Element 


OescHotlon:  Defines  a  uniform  static  pressure  load  applied  to  two-dimensional  elements.  Only 
UURIJ3 , ""SHEER ,  TRIM,  or  TRIAX6  elements  may  have  a  pressure  load  applied  to  them  via  this  card. 


Format  and  Example: 


Alternate  Form: 


*  P  Pressure  value  (Real) 

EIO  ) 

EID1>  Element  Identification  number  (Integer  >  0;  EID1  <  EID2) 

EI02J 

Remarks:  1.  EIO  must  be  0  or  blank  for  omitted  entries. 

2.  Load  sets  must  be  selected  In  the  Case  Control  Oeck  (L0AD»SIO)  to  be  used  oy 
MSC/NASTRAN. 

3.  At  least  one  positive  EID  must  be  present  on  each  PL0AD2  card. 

4.  If  the  alternate  form  Is  used,  all  elements  EID1  through  EID2  must  be  two- 

dimensional. 

5.  The  direction  of  the  pressure  is  computed  according  to  the  right-hand  rule  using 

the  grid  point  seouence  specified  on  the  element  card.  Refer  to  the  PL0AD  card. 

6.  All  elements  referenced  must  exist. 

7.  Continuation  cards  are  not  allowed. 
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Input  Data  Card  NLPARM  Parameters  for  Nonlinear  Analysis  Control 
Description:  Defines  a  set  of  parameters  for  nonlinear  analysis  iteration  strategy 
Format  and  Example: 


1  2  3  4  5  6  7  8  9  10 


NLP ARM 

ID 

INC 

OT 

KM6TH0D 

KSTEP 

MAX  I TER 

C0NV 

INT3UT 

NLPARM 

15 

5 

PV 

+NL1 

EPSU 

EPSP 

EPSW 

MAXOIV 

MAXQN 

MAXLS 

FSTRESS 

*NL1 

1.0E-3 

1.0-4 

1.0-5 

1 

20 

8 

.02 

Contents 

Identification  Number  (Integer  >  D) 

Number  of  Increments  (Integer  >  0) 

Incremental  time  Interval  for  creep  analysis  (Real  _>  0)  (Oefault  ■  0.;  no  creep) 

Method  for  controlling  tangent  stiffness  updates  (BCO  •  AUT0, ITER, SEMI .LSQN.AUT0QN, 
or  SEMIQN)  (No  default).  See  Remark  S. 

Number  of  Iterations  before  a  matrix  update.  Used  only  for  ITER  and  LSQN 
methods.  (Integer^O  or  blank)  (Default  *  I) 

Limit  on  total  Iterations  for  each  load  increment  (Integer  k  0)  (Default  »  20) 

Flags  to  select  convergence  tests  (BCO  •  U,  P,  W,  or  any  combination  of  the 
letters)  (Oefault  *  W).  See  Remarks  3  and  4. 

Intermediate  output  flag  (BCO  «  TES  or  N0)  (Default  •  YES) .  See  Remark  5. 

Error  tolerance  limits  for  determining  convergence  (Real  >  0.0)  (Defaults  ■  1.0E-2, 
1.0E-3,  1.0E-4,  respectively) 

Limit  on  diverging  Iterations  (Integer  >  0)  (Oefault  •  l) 

Maximum  number  of  Quasi -Newton  correction  vectors  to  be  saved  on  the  data  base 
(Integer  0)  (Oefault  ■  10) 

Maximum  number  of  linear  search  operations  per  vector  Iteration  (Integer  >  0) 
(Oefault  -  3) 

Fraction  of  effective  stress  (7)  used  to  determine  the  subincrement  size  in  the 
material  routines  (0.  <  Real  <  1.0)  (Default  •  0.15).  See  Remark  9. 


Field 

10 

INC 

OT 

KMETH0P 

KSTEP 

MAX  ITER 
C0NV 

INT0UT 

EPSU.EPSP. 

EPSM 

MAXOIV 

MAXQN 

MAXLS 

FSTRESS 
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ML P ARM  (Cont.) 


The  MLPARM  8ulk  Data  cart  Is  requested  by  the  Case  Control  cart  NLP  ARM  *  ID.  Each 
solution  subcase  requires  a  loading  condition  and  an  MLPARM  request. 

In  case  of  static  analysis  (QT  •  blank),  CMC  Is  the  number  of  equal  subdivisions  of 
the  net  load  change  defined  for  the  SUBCASE.  Applied  loads,  gravity  loads, 
temperature  sets,  enforced  displacements,  etc.,  define  the  new  loading  condition. 
The  differences  from  the  previous  case  are  divided  by  INC  to  define  the  incremental 
values.  In  case  of  creep  analysis  (DT  >  0.),  INC  1$  the  number  of  time  step 
Increments. 

The  test  flags  (U.P.W)  and  the  error  limits  (EPSU,  etc.)  define  the  convergence 
criteria  (U  *  Displacement  error  test,  P  •  Load  equilibrium  error  test,  W  *  work 
error  test).  If  the  Internally  calculated  error  fractions  are  less  than  the 

requested  limits,  the  Iteration  stoos,  the  results  are  processed,  and  the  program 
continues  to  the  next  load  increment  or  subcase.  All  requested  error  criteria  must 
be  satisfied  for  a  "convergent  solution*. 

If  convergence  does  not  occur  for  a  particular  load  Increment,  the  number  of 
Iterations  Is  limited  to  MAXITER.  If  MAXITER  is  a  positive  number,  the  results 

will  be  processed  and  the  program  will  continue  with  the  next  load  increment.  If 

MAXITER  Is  a  negative  value,  the  solution  will  terminate  if  convergence  Is  not 

achieved. 

The  basic  nonlinear  solution  methods  are  the  “QFGS"  Quasi-Newton  and  modified 
"Newton-Raphson"  processes  whereby  the  out-of -balance  nonlinear  loads  are  measured 
and  recycled  to  solve  for  Improved  displacements.  The  convergence  of  this 
iteration  process  is  generally  Improved  by  efther  performing  'line  searches*  or  by 
frequent  updates  to  the  "tangent  stiffness  matrix*  to  account  for  current  plastic 
and  geometric  effects.  However,  the  matrix  updates  may  be  relatively  costly  and 
may  be  detrimental  if  the  estimates  are  far  from  the  correct  answer.  Thus,  the 
user  Is  given  the  following  (METHOD  options  to  provide  control  over  a  large  range 
of  problem  types. 

a.  If  the  AUTO  option  is  selected,  the  program  will  automatically  select  the  most 
efficient  strategy  based  upon  convergence  rates.  At  each  step  the  number  of 
iterations  required  to  converge  Is  estimated.  If  this  exceeds  the  MAXITER 
limit,  or  If  the  CPU  time  exceeds  the  time  to  perform  a  stiffness  matrix 
update,  the  module  will  perform  a  stiffness  update.  If  divergence  occurs  for 
MAXOIV  successive  iterations,  a  stiffness  update  is  also  performed. 

b.  If  the  SEMI  option  Is  selected,  for  each  load  increment  the  program  will;  (I) 
perform  a  single  displacement  iteration  based  upon  the  new  load;  (2)  update 
the  stiffness  matrix  for  the  estimated  geometry  and  material  properties;  and 
(3)  resume  the  normal  AUTD  Iteration  method. 

c.  If  the  ITER  option  Is  selected,  the  stiffness  matrix  will  be  updated  only 
after  each  (STEP  Iteration.  The  count  Is  reset  to  zero  for  each  new  matrix  or 
load  step. 

d.  If  one  of  the  AUT0QN,  SEMIQN,  or  the  ISQN  options  Is  selected  the  “Quasi  - 
Newton*  and  “Line  Search*  procedures  are  used.  Basically  the  line  search  (LS) 
will  scale  the  displacement  increment  vector  to  minimize  the  error  at  each 
Iteration.  If  the  line  search  provides  a  substantial  improvement,  the  Quasi- 
Newton  (QN)  operations  use  the  local  tangent  of  the  error  to  correct  the 
tangent  matrix  solution  for  subsequent  iterations.  Otherwise  the  operations 
behave  much  like  the  AUTO,  SEMI,  and  ITER  options. 


(Continued) 


NLP ARM  (Cont.) 


6.  Output  requests  for  ELF0RCE  and  STRESS  mad*  in  Case  Control  nill  be  processed  if 

1NT0UT  »  YES  (the  default  value).  If  INT0UT  •  NO,  these  output  requests  are 

processed  for  only  the  last  load  factor  In  a  subcase.  Requests  for  OISP  printout 

and  structure  plots  are  also  processed  during  the  nonlinear  Iteration  if  PARAM, 

NLOISP ,  1  Is  provided  In  the  Bulk  Data  deck. 

7.  If  a  diverging  solution  is  encountered  for  MAXOIV  successful  iterations,  the  AUTO, 

AUT0QN,  SEMI,  and  SEMIQN  options  Mill  perform  a  tangent  stiffness  update.  If, 

after  the  stiffness  update,  the  solution  Is  still  diverging,  the  program 
terminates.  For  ITER  and  ISQK  options,  diverging  cases  cause  immediate 
termination. 

8.  The  unit  of  the  incremental  time  interval  under  OT  field  must  be  consistent  with 
the  unit  used  on  the  CREEP  card  to  define  the  creep  characteristics.  Total  creep 
time  for  Che  subcase  may  be  found  from  OT  multiplied  by  the  value  in  the  INC. 

9.  The  number  of  subincrements  In  the  material  routines  (elastoplastlc  and  creep)  are 
determined  such  that  the  subincrement  size  is  approximately  FSTRESS  *  7. 

FSTRESS  is  also  used  to  establish  a  tolerance  for  elastoplastlc  material  in  the 
material  routine,  i.e.. 


4  in  yield  function  <  FSTRESS  •  yield  stress 


If  the  limit  Is  exceeded  at  the  converging  state,  the  program  Mill  exit  elth  a 
fatal  error  message.  Otherwise  the  stress  state  is  adjusted  to  the  current  yield 
surface,  resulting  in  4  ■  0. 


m 


''i  4  •  *  .  **•  ■  m  *’»  **#  •*  #  •**  f  " 


Input  Oat*  Card  SPCS  Constraint  Sat  Sanarator,  Standard  Form 

m 

Oascrlptlon:  Sancratas  SPC1  cards  that  constrain  a  subset  of  grid  points  In  a  grid  point  field 
Format  and  Example; 


.  * *  -  . 


Input  Data  Card  OLgAO  Dynamic  Load  Combination  (Superposition) 

Description:  Defines  a  dynamic  loading  condition  for  frequency  response  or  transient  response 
problems  "as  a  linear  combination  of  load  sets  defined  via  RL3AD1  or  RL8A0Z  cards  (for  frequency 
response)  or  TL3AD1  0r  TL3ADZ  cards  (for  transient  response) 

Format  and  Examples : 


0L3AD 

17 

1.0 

2. 

Contents 

Load  set  identification  number  (Integer  >  0) 

Scale  factor  (Real) 

Scale  Factors  (Real) 

Load  set  identflcation  numbers  defined  via  card  types  enumerated  above 
(Integer  >  0) 


Remarks:  1.  The  load  vector  being  defined  by  this  card  is  given  by 


W  •  S  I  SjIM  . 
1  1  1 


The  Li  must  be  unique. 

SID  must  be  unique  from  all  Li. 

Nonlinear  transient  loads  may  not  be  included;  they  are  selected  separately  in  the 
Case  Control  Deck. 

Linear  load  sets  nust  be  selected  in  the  Case  Control  Deck  (DL3A0»SI0)  to  be  used 
by  MSC/NASTRAN. 

A  OLJAO  card  may  not  reference  a  set  identification  number  defined  by  another  DL3A0 
card. 

TL3AQ1  and  TL3A02  loads  may  be  combined  only  through  the  use  of  the  OUJAO  card. 
RL3A01  and  R19A02  loads  may  be  combined  only  throu^i  the  use  of  the  OL3AO  card. 

SID  must  be  unique  for  all  TL3AD1 ,  TL0AO2,  RL3A01  and  RL3AD2  cards. 


Input  Oata  Card  TlgAPl  Transient  Response  Dynamic  Load,  Form  I 

Description:  Oeflnes  a  time  dependent  dynamic  load  or  enforced  motion  of  tfte  form 


•  (A  F(t  -  t)l 


for  use  In  a  transient  response  proDlt 


Format  and  Example: 


TUA01 


TL3AOI 


Contents 

Set  Identification  number  (Integer  >  0) 


Identification  number  of  OAREA  card  set  or  a  static  load  card  set  or  a  thermal  load 
set  (In  heat  transfer  analysis)  which  defines  A  (Integer  >  0) 


Identification  number  of  OEIAY  card  set  which  defines  r  (Integer  _>  0) 


Defines  the  nature  of  the  dynamic  excitation  (Integer  0,  1,  2,  3  or  blank).  (See 
Remark  2.) 


Identification  number  of  TA8LE01  card  which  gives  F(t  -  t)  (Integer  >  0) 


Remarks:  1.  If  M  Is  zero,  t  will  be  zero. 


2.  The  nature  of  the  dynamic  excitation  Is  defined  In  field  5  In  accordance  with  the 
following  table. 


Integer 

Excitation  Function 

0  or  blank 

Force  or  Moment 

l 

Enforced  Displacement 

2 

Enforced  Velocity 

3 

Enforced  Acceleration 

See  Section  2.7  of  the  MSC/NASTRAR  Application  Manual  regarding  the  use  of  the 
enforced  motion  options.  Note  that  "large  masses*  must  be  used  for  enforced 
motion.  For  heat  transfer  problems.  Field  S  must  be  blank. 


(Continued) 
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TUJA01  (Cont.) 


3.  Dynamic  load  sets  mist  be  selected  In  the  Case  Control  Qecic  (DL0AO*SIO)  to  be  used 
by  MSC/NASTRAN. 

4.  TL0AO 1  loads  may  be  combined  with  TL3A02  loads  only  by  soeclfl cation  on  a  0L3A0 
card.  That  Is,  the  SID  on  a  T13A01  card  may  not  be  the  same  as  that  on  a  TL0AO2 
card. 

5.  SID  must  be  unique  Tor  all  TL0AO1,  TL0AO2,  RL0AOI  and  RL0AD2  cards. 

6.  Field  3  may  reference  sets  containing  QHBOY,  qBOYl,  QB0Y2,  QVECT,  and  QV0L  cards 

when  using  the  heat  transfer  option. 

7.  If  the  heat  transfer  option  Is  used,  the  referenced  QVECT  data  card  may  also 

contain  references  to  functions  of  time,  and  therefore  A  may  be  a  function  of  time. 

8.  Fourier  analysis  will  be  used  If  this  Is  selected  In  an  aeroelastlc  response 

problem. 


Input  Data  Card  LSEQ  Static  Load  Set  Definition 


Description:  Defines  a  seouence  of  static  load  sets  to  be  applied  to  the  structural  model.  The 
load  sets  may  be  referenced  by  dynamic  load  cards. 


Format  and  Example: 


Field  Contents 

SID  Set  Identification  of  the  set  of  LSEQ  cards  (Integer  >  0) 

OAREA  The  0AREA  set  identification  assigned  to  this  static  load  vector  (Integer  >  0) 

LID  Load  set  Identification  number  of  a  set  of  static  load  cards  (Any  card  that  may  be 

referenced  by  the  L0AO  Case  Control  card)  (Integer  >  0  or  blank) 

TID  Set  identification  of  a  thermal  load  set  (Any  card  that  may  be  referenced  by  the 

TEMP(LWO)  Case  Control  card)  (Integer  >  0  or  blank) 


Remarks:  I.  The  above  cards  will  not  be  used  unless  selected  in  the  Case  Control  Deck  with  a 
L0ADSET  card. 

2.  This  card  Is  available  only  In  suoerel ements  and  In  Solution  Sequences  26,  27,  30 
and  3 i. 

3.  The  number  of  static  load  vectors  created  for  each  suoerel ement  is  the  number  of 
unlaue  OAREA  10s  on  all  LSEQ  cards  In  the  bulk  data. 

4.  The  OAREA  Identification  assigned  to  the  static  load  vectors  may  be  referenced  by 
RL3A01.  RL0AO2,  TL3A01  and  TL3A02  cards. 

5.  Element  data  recovery  for  thermal  loads  Is  not  currently  Imolemented  in  dynamics. 

6.  OAREA  set  Identification  numbers  should  be  unique  with  respect  to  all  static  'oao 
set  Identification  numbers. 
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UNCLASSIFIED 


F/G  19/4 


Input  Oata  Card  TABLED I  Oynanric  Load  Tabular  Function,  Form  1 


Description:  Defines  a  tabular  function  for  use  In  generating  frequency-dependent  and  time- 

dependent  dynamic  loads 


Format  and  Example: 


1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

TABLED! 

10 

x 

lX 

xx 

X 

XX  1 

TABLE01 

32 

| 

♦ABC  ! 

*1 

n 

*2 

*2 

*3 

*3 

*4  I 

n 

1 

♦8C 

-3.0 

S.9 

2.0 

5.6 

3.0 

5.6 

ENOT 

-  etc.- 


Fleld  Contents 

ID  Table  identification  number  (Integer  >  0) 


*1.*1 


Tabular  entries  (Real) 


Remarks:  l.  The  x1  must  be  In  either  ascending  or  descending  order  but  not  both. 

2.  Jumps  between  two  points  (xf  >  xul)  are  allowed,  but  not  at  the  end  points. 

3.  At  least  two  entries  mist  be  present. 


♦.  Any  x-y  entry  nay  be  Ignored  by  placing  the  BCD  string  "SKIP*  In  either  of  the  two 
fields  used  for  that  entry. 

5.  The  end  of  the  table  Is  Indicated  by  the  existence  of  the  BCD  string  "ENOT"  in 
either  of  the  two  fields  following  the  last  entry.  An  error  Is  detected  if  any 
continuation  cards  follow  the  card  containing  the  end-of -table  flag  "ENOT*. 

S.  Each  TABLE01  mnemonic  Infers  the  use  of  a  specific  algorithm.  For  TA8LED1  type 
tables,  this  algorithm  is 


▼  ■  yt(x) 


where  X  Is  Input  to  the  table  and  Y  is  returned.  The  table  look-up  yT(xj,  x  *  X, 
Is  performed  using  linear  Interpolation  within  the  table  and  linear  extrapolation 
outside  the  table  using  the  last  two  end  points  at  tne  approorlate  table  end.  At 
jum>  points  the  average  y-(x)  is  used.  There  are  no  error  returns  from  this  table 
look-up  procedure. 

7.  Linear  extrapolation  Is  not  used  for  Fourier  Transform  methods.  The  function  1$ 
zero  outside  the  range. 


Input  Oats  Card  TSTEPNL 


Tina  Step  and  Nonlinear  Solution  Controls 


Description:  Provides  parametric  controls  and  data  for  Nonlinear  Transient  Analysis 
format  and  Example: 

1  2  3  4  5  6  7  8 


10 


TSTEP NL 

ID 

NOT 

DT 

NO 

METHOD 

KSTEP 

MAXI  TER 

CDNV 

tstepnl 

250 

40 

.0015 

5 

AUTO 

5 

P 

+NL1 

EPSU 

EPS? 

EPSU 

MAXOIV 

MAXQN 

MAXLS 

PSTRESS 

♦NL1 

.01 

1 

20 

6 

.02 

field  Contents 

10  Identification  number  (Integer  >  0) 

NOT  Number  of  time  steps  of  value  0T(1)  (Integer  _>  2) 

OT  Time  Increment  (Peal  >  0.0) 

NO  Skip  factor  for  output  (every  NQth  step  will  be  saved  for  output)  (Integer  >  0) 

METHOD  Method  for  controlling  Iterations  and  tangent  matrix  updates  (BCD  •  AUTO  or  TSTEP ) 

(No  default) 

A  STEP  Number  of  steps  before  a  matrix  update  (Integer  ^  0)  (Oefault  ■  0) 

MAXITER  Limit  on  total  Iterations  for  each  time  increment  (Integer  >  0)  (Default  •  2) 

CONV  flags  to  select  convergence  tests  (BCD  »  U,  P,  W,  or  any  combination  of  the 

letters  or  blank) 

EPSU  ,EPSP  .EPSU  Error  tolerance  limits  for  determining  convergence  (Real  >  0.0) 

(Oefaults  •  l.QE-2,  l.OE-3,  1.0E-4,  respectively) 

NAXDIV  Limit  on  the  number  of  diverging  Iterations  at  each  time  and  on  diverging  time 

steos  (Integer  >  0)  (Default  ■  2) 

maxqn  Maximum  number  of  Quasi -Newton  correction  vectors  to  be  saved  on  the  data  base 

(Integer  >  0)  (Oefault  ■  20) 

NAXLS  Maximum  number  of  Line  Search  Iterations  for  LSQN  option  (Integer  >  0) 

(Oefault  -2)  ~ 

PSTRESS  fraction  of  effective  stress  f«)  used  to  determine  the  subincrement  size  in  the 

material  routines  (0.  <  Reel  <  1.0)  (Default  *  0.15).  See  Remark  10. 
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TSTEPNL  (Cont.) 


Remarks:  1.  Tti«  TSTEPNL  Sulk  0*ta  card  Is  selected  by  the  Cast  Control  card  TSTEPNL  ■  ID.  Each 
suOcast  (residual  superelement  soludons  only)  requires  a  TSTEPNL  card  and  either 
applied  loads  via  TL0AO1  data  or  initial  values  froa  a  previous  subcase.  Multiple 
subcases  are  assuned  to  occur  sequentially  In  tine  with  the  Initial  values  of  time 
and  displacement  conditions  of  each  subcase  defined  by  the  end  conditions  of  the 
previous  subcase. 

2.  OT,  NOT,  N0  are  the  time  Increments,  at,  the  total  number  of  steps  Nt,  and  the  out* 
put  Interval,  ng,  respectively.  The  value  of  time,  Tn,  at  each  Increment  n  will  Be 


tfl  -  tg  *  nAt  n  •  0,  I.....N, 


where  tg  Is  the  Initial  value  of  time  which  corresponds  to  the  last  step  from  the 
previous  subcase. 

3.  For  printing  and  plotting  the  solution,  data  will  be  output  at  steps  n  ■  0,  ng, 

2ng,...N.  The  Case  Control  card  JTIME  may  also  be  used  to  control  the  output 

points. 

4.  The  basic  solution  method  Is  the  "Newmarfc  Beta*  transient  numerical  method  along 

with  the  "Newton-Raphson”  method  for  controlling  nonlinear  effects.  Equilibrium 

solutions  with  respect  to  Inertial,  damping,  elastic,  and  nonlinear  loads  are 
obtained  at  each  time  step  much  the  same  as  In  the  statics  case. 

The  KMETHUO  options  provide  the  user  with  three  optional  nonlinear  Iteration 

methods.  These  are: 

AUTJ  provides  for  automatic  stiffness  matrix  updates  to  improve  poor 
convergence.  The  KSTEP  value  Is  Ignored. 

TSTEP  provides  for  a  matrix  update  at  every  KSTEPth  Increment  of  time 
(typically  a  large  number). 

5.  The  stiffness  matrix  Is  always  updated  for  a  new  subcase  or  restart.  Only  one 

update  Is  allowed  for  each  time  step  and  the  Iteration  process  returns  to  the 
previous  time  step  after  an  update. 

6.  MAX  I TER  defines  the  number  of  static  nonlinear  load  Iterations  (corrections)  to  be 
performed  at  each  time  step.  At  least  two  iterations  are  necessary  to  determine 
the  convergence  of  the  displacements  (see  below). 

7.  The  test  flags  (U,P,K)  and  tne  error  limits  (EPSU,  etc.)  define  the  convergence 

criteria  (U  •  Displacement  error  test,  P  *  Load  equilibrium  error  test,  w  »  work 
error  test).  If  the  Internally  calculated  error  fractions  are  less  than  the 

requested  limits,  the  Iteration  stops,  the  results  are  processed,  and  the  program 
continues  to  the  next  time  step  or  subcase. 


TSTEPNt.  (Cont.) 


8.  If  no  Costs  sro  requested  or  no  convergence  occurs,  the  number  of  Iterations  is 

Halted  to  MAXITER.  The  results  will  be  processed  and  the  program  will  continue 

Mlth  the  next  tine  step. 

9.  MAXQN  provides  the  control  over  diverging  solutions.  At  each  time  step,  the 

solution  Is  allotted  to  diverge  until,  at  the  MAXOIV  iteration,  a  matrix  update  is 

requested.  If  the  average  convergence  parameter  indicates  divergence  for  HAXOIV 

successive  time  steps,  the  transient  solution  stops  and  exits  for  data  recovery. 

10.  The  number  of  sublncrements  in  the  material  routines  (elastoplastlc  and  creep)  are 
determined  such  that  the  subtnerement  size  Is  approximately  FSTRESS  *  7. 

F STRESS  Is  also  used  to  establish  a  tolerance  for  elastoplastlc  material  in  the 
material  routine,  l.e.. 


4  in  yield  function  <  FSTRESS  *  yield  stress 


If  the  Unit  Is  exceeded  at  the  converging  state,  the  program  will  EXIT  with  a 
fatal  error  massage.  Otherwise  the  stress  state  is  adjusted  to  the  current  yield 
surface,  resulting  In  4  ■  0. 
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